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[i] Ozone changes and associated climate impacts in the Coupled Model Intercomparison 
Project Phase 5 (CMIP5) simulations are analyzed over the historical (1960-2005) and 
future (2006-2100) period under four Representative Concentration Pathways (RCP). In 
contrast to CMIP3, where half of the models prescribed constant stratospheric ozone, 

CMIP5 models all consider past ozone depletion and future ozone recovery. Multimodel 
mean climatologies and long-term changes in total and tropospheric column ozone 
calculated from CMIP5 models with either interactive or prescribed ozone are in 
reasonable agreement with observations. However, some large deviations from 
observations exist for individual models with interactive chemistry, and these models are 
excluded in the projections. Stratospheric ozone projections forced with a single halogen, 
but four greenhouse gas (GHG) scenarios show largest differences in the northern 
midlatitudes and in the Arctic in spring (~20 and 40 Dobson units (DU) by 2 1 00, 
respectively). By 2050, these differences are much smaller and negligible over Antarctica 
in austral spring. Differences in future tropospheric column ozone are mainly caused by 
differences in methane concentrations and stratospheric input, leading to -10DU increases 
compared to 2000 in RCP 8.5. Large variations in stratospheric ozone particularly in 
CMIP5 models with interactive chemistry drive correspondingly large variations in lower 
stratospheric temperature trends. The results also illustrate that future Southern Hemisphere 
summertime circulation changes are controlled by both the ozone recovery rate and the rate 
of GHG increases, emphasizing the importance of simulating and taking into account 
ozone forcings when examining future climate projections. 

Citation; Eyring, V., et al. (2013), Long-term ozone changes and associated climate impacts in CMIP5 simulations, 
J. Geophys. Res. Atmos., 118 , 5029-5060, doi:10.1002/jgrd.50316. 


1. Introduction 

[ 2 ] To improve understanding of climate and to provide 
estimates of fuUire climate change, the Working Group 
on Coupled Modelling (WGCM) of the World Climate 
Research Programme (WCRP) has defined coordinated 
climate model experiments, fonning phase 5 of the Coupled 
Model Intercomparison Project (CMIP5) [Taylor et al., 
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2012], Since the previous round (CMIP3), the participating 
models have been further developed, new models have 
joined, and new emission scenarios, so-called Representa- 
tive Concentration Pathways (RCPs), have been established 
[Meinshausen et al., 2011; Moss et al., 2010]. The focus of 
this study is to examine the stratospheric and tropospheric 
ozone evolution in the CM1P5 historical (1960-2005) and 
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future (2006-2100) simulations under varying greenhouse 
gas (GHG) and halogen concentrations and to highlight 
some related climate impacts. 

[ 3 ] Stratospheric ozone depletion has been found to play a 
prominent role in recent Southern Hemisphere (SH) climate 
change, with impacts on the tropospheric circulation and 
surface climate in austral summer [Kang et al., 2011; Son 
et al., 2010; Thompson et al., 2011], For future projections, 
the interplay between stratospheric ozone recovery from a 
decrease of ozone-depleting substances (ODSs) [ WMO , 
2011] and increasing GHGs could in addition have impor- 
tant long-tenn implications for CCF uptake in the Southern 
Ocean [Le Quere et al., 2007; Lenton et al., 2009] and 
changes to Antarctic ice sheets [Turner et al., 2009a], 
Recent studies have also suggested that changes in strato- 
spheric ozone could have potential impacts on SH sea ice 
[Bitz and Polvani, 2012; Sigmond and Fyfe, 2010; Smith 
et al, 2012; Turner et al., 2009b]. To date, stratospheric 
ozone projections have been produced with Chemistry- 
Climate Models (CCM) [WMO, 2007; 2011], where sea 
surface temperatures (SSTs) and sea-ice concentrations (SICs) 
are mostly prescribed from observations or climate simula- 
tions. Similar to CM1P, coordinated community-wide simula- 
tions for stratospheric ozone have been organized by the 
Chemistry-Climate Model Validation (CCMVal) Activity of 
WCRP’s Stratospheric Processes and their Role in Climate 
(SPARC) project [Austin et al., 2010; Eyring et al., 2010b; 
Eyring et al., 2006; Eyring et al., 2007; SPARC-CCMVal, 
2010], In these stratospheric-resolving CCMs, chemistry is 
interactively coupled to the dynamics and radiation, allowing 
for chemistry-climate feedbacks. Increasingly, the chemistry 
and dynamics of the stratosphere and troposphere are 
now being modeled as a single entity. The 2010 WMO/UNEP 
Scientific Assessment of Ozone Depletion [WMO, 2011] 
featured several stratospheric models that included tropo- 
spheric chemistry, and one model with a coupled ocean. 
It was one of the main recommendations of the SPARC- 
CCMVal [2010] report that stratosphere-resolving CCMs 
should continue to evolve towards more comprehensive, 
self-consistent stratosphere-troposphere CCMs, including a 
coupled ocean. 

[ 4 ] In CMIP3, which supported the Intergovernmental 
Panel on Climate Change (IPCC) Fourth Assessment 
Report (ART) [IPCC, 2007] with climate simulations, all of 
the ~20 models examined for example by Son et al. [2010] 
used prescribed ozone. However, ozone in the CMIP3 simu- 
lations is not documented in detail, since it was not archived 
from the model simulations, and only some rough infonnation 
exists. Around half of the CMIP3 models prescribed a time- 
varying stratospheric ozone tune series, for example scaled 
by the evolution of equivalent effective stratospheric 
chlorine (EESC). The other half prescribed a constant ozone 
climatology, thus not accounting for climate effects induced 
by stratospheric ozone depletion in the past and projected 
recovery in the future. This led to biases in the tropospheric 
general circulation of the SH compared to observations 
[Cordero and Forster, 2006], The missing Antarctic ozone 
depletion in late austral spring resulted in weaker poleward 
displacement of the tropospheric midlatitude jet and smaller 
expansion of the SH Hadley cell in summer \Son et al., 
2009; Son et al., 2010] than in models with prescribed time- 
varying ozone. For the future projections, the CM1P3 models 


without ozone recovery did not simulate a reversal of the 
zonal mean wind trends that was projected over the first half 
of the 21st century by both CCMs with interactive strato- 
spheric chemistry [Perlwitz et al., 2008; Son et al., 2008; 
Son et al., 2010; SPARC-CCMVal, 2010; WMO, 2011] and 
CM1P3 models that prescribed ozone recovery [Son et al., 
2010], In addition, the incomplete representation of 
stratospheric ozone leads to biases in temperature trends in 
the upper troposphere/lower stratosphere (UTLS) region 
[Forster et al., 2007; Lanzante, 2007; Polvani and Solomon, 
2012; Santer et al., 2012; Solomon et al., 2012; Young 
et al., 2013a], 

[ 5 ] Similarly, chemistry-climate coupling and the consid- 
eration of time-varying ozone are also important in the 
troposphere, where ozone is of concern as a pollutant and 
GHG. Tropospheric ozone concentrations depend on the 
magnitude and spatial distribution of the emissions of its 
precursors, nitrogen oxides (NO x ), carbon monoxide (CO), 
and volatile organic compounds (VOCs), as well as on 
climate. Several studies have investigated the impact of 
emissions and climate changes on future ozone for a range 
of scenarios. In a warming climate, tropospheric ozone loss 
can increase due to higher absolute humidity [Johnson 
et al., 1999; Zeng et al., 2008]. On the other hand, climate 
change may increase tropospheric ozone, for example, 
through an increased influx of stratospheric ozone [Collins 
et al., 2003; Hegglin and Shepherd, 2009; Sudo et al., 
2003] or through increases in lightning NO x emissions 
[Price and Rind, 1994; Schumann and Huntrieser, 2007], 
although static [Stevenson et al., 2006] or decreasing 
[Jacobson and Streets, 2009] emissions have also been 
projected. Climate-dependent emissions of biogenic VOCs 
[Arneth et al., 2010; Guenther, 2007; Guenther et al., 
2006] are another potential climate feedback, although 
whether they will increase or decrease in the future is unclear 
[Arneth et al., 2010; Hauglustaine et al., 2005; Lathiere 
et al., 2005; Sanderson et al., 2003; Young et al., 2009]. 
While a positive relationship between anthropogenic emis- 
sions and tropospheric ozone has been shown [Stevenson 
et al., 2006; Young et al., 2013b], the net impact of climate 
and emission changes is uncertain and is likely to vary sig- 
nificantly by region, altitude, and season [Isaksen et al., 
2009; Jacob and Winner, 2009; Stevenson et al., 2006]. 

[6] As a response to the inconsistencies in ozone 
representation in the CMIP3 models, the CMIP5 protocol 
recommended a time-varying ozone database for models 
without interactive chemistry. This continuous tropospheric 
and stratospheric vertically resolved ozone time series from 
1850 to 2099 was developed as a joint effort of the Interna- 
tional Global Atmospheric Chemistry (IGAC) and SPARC 
projects, hereafter referred to as the IGAC/SPARC ozone 
database [Cionni et al., 2011], The historical stratospheric 
zonal mean data set is based on the Randel and Wu [2007] 
observational database from 1979 to 2009 and extrapolated 
backward to 1850 with the help of a regression representing 
EESC and the 1 1 year solar cycle variability. The data set 
continues beyond 2009 using the multimodel mean 
projections from 13 CCMs that performed a simulation 
under the SRES (Special Report on Emission Scenarios) 
A1B GHG scenario [IPCC, 2000] and the Al adjusted 
halogen scenario [WMO, 2007] in the second round of 
CCMVal (CCMVal-2). The stratospheric zonal mean ozone 
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time series is merged with a three-dimensional (3-D) tropo- 
spheric data set extracted from simulations of the past by two 
CCMs (CAM3.5 [ Lamarque et al., 2012] and Goddard 
Institute for Space Studies (GISS)-PUCCINI [Shindell et al., 
2006]) and of the future by one CCM (CAM3.5). The future 
tropospheric ozone time series continues the historical 
CAM3.5 simulation until 2099 following the four different 
RCPs. The majority, but not all, of the CM1P5 models with- 
out interactive chemistry prescribed ozone according to this 
or a modified version of the IGAC/SPARC ozone database 
(see Table 1 and details in the Appendix in section A2.1). 
However, a subset of nine of the CMIP5 models also includes 
interactive chemistry, spanning the surface through the 
stratosphere and above. 

[ 7 ] This study thus represents the first opportunity to 
undertake a multimodel assessment of ozone and associated 
climate impacts in fully coupled transient chemistry-climate 
simulations and in the new set of CMIP5 models. We 
first evaluate the historical CMIP5 ozone simulations with 
observations, then assess the ozone projections, and finally 
highlight prominent features of associated climate impacts. 
In addition, we present a detailed documentation of ozone 
in the CM1P5 simulations both for the multimodel mean 
(or different classes of multimodel means, see below) as 
well as for individual models. We therefore provide a 
description of the chemistry schemes or the prescribed 
ozone database for each CMIP5 model in the Appendix 
and include individual models in the evaluation of ozone. 
The analysis presented here follows on from a previous 
study [. Eyring et al., 2010a], which examined differences 
in future stratospheric ozone for six different GHG scenarios 
but only with a very limited set of models. In particular, 
there was only one model included that simulated ozone 
under the RCP scenarios. Tropospheric ozone is also briefly 
discussed to show that it contributes substantially to the 
differences in total column ozone among the four RCPs in 
the 21st century, complementing Young et al. [2013b] who 
presented the evolution of tropospheric ozone with the 
historical and RCP scenarios in a set of (mostly) atmosphere- 
only models that participated in the Atmospheric Chemistry 
Climate Model Intercomparison Project (ACCMIP). 

[s] The paper is organized as follows. Section 2 describes 
the CMIP5 models and simulations, with a focus on how 
ozone chemistry was treated (details are presented in the 
Appendix). In section 3, ozone in the historical simulations 
is evaluated against observations, and projections of ozone 
in the 21st century are provided for the four RCP scenarios, 
using only the group of CMIP models with interactive or 
semi-offline chemistry (see section 2 for details). In 
section 4, climate impacts associated with stratospheric 
ozone changes are discussed for the past and future. The 
focus of the analysis is on different groups of multimodel 
mean analysis (see section 2), with the performance of 
individual CMIP5 models being discussed as part of the 
model evaluation. We close with a summary in section 5. 

2. Models and Simulations 

[ 9 ] In this study, the CMIP5 historical (1850-2005) and 
future simulations under the RCPs (2006-2100) [Taylor 
et al., 2012] are analyzed for a total of 46 different model 
versions from 22 different modeling centers; see Table 1. 


The model simulations are driven by concentrations of 
GHGs and also include changes in aerosols, ozone, land 
cover, and natural forcings. In order not to bias the 
multimodel mean towards the available number of ensem- 
bles from each model version, only one ensemble member 
is considered. It should be noted that the effective number 
of independent models in CMIP5 is smaller [Annan and 
Hargreaves, 2011; Jun et al., 2008; Knutti et al., 2013; 
Masson and Knutti, 2011; Pennell and Reichler, 2011]. 

[ 10 ] Since future ozone projections can only be derived 
from models with interactive chemistry, the models are 
grouped into two different categories: the first group 
includes CMIP5 models with interactive or semi-offline 
chemistry (CHEM, see section 2.1), and the second group 
includes models with prescribed ozone (NOCHEM, see 
section 2.2). The multimodel means are shown for these 
two groups as well as for the multimodel mean calculated 
from all CMIP5 simulations considered. Individual models 
are shown in selected figures in the main part of the paper 
and in addition in the Appendix. 

[ 11 ] The main difference between the two groups is for 
stratospheric ozone, where the CHEM models were run 
under the emission specifications for each of the RCP 
scenarios, but the NOCHEM simulations prescribed ozone 
from simulations that were driven only with a single GHG 
scenario. The prescribed data set was either based on the 
original or a slightly modified version of the IGAC/SPARC 
ozone database or an alternate data set following just a 
single GHG scenario. Therefore, future stratospheric ozone 
concentrations are identical in all four RCPs in the 
NOCHEM group (see section 2.2). The ODSs are prescribed 
according to the Al scenario of WMO [2003] in all RCPs 
[Meinshausen et al., 2011], This scenario does not include 
the earlier phase out of hydrochlorofluorocarbons that 
was agreed to by the Parties to the Montreal Protocol in 
2007. It is therefore slightly different from the adjusted 
halogen scenario of WMO [2007] that was used in the 
CCMVal-2 simulations [SPARC-GCMVal, 2010; WMO, 
2011] on which the IGAC/SPARC ozone data set is based. 
Ozone in the troposphere however follows ozone precursor 
emissions under the four RCPs in all CMIP5 simulations 
(CHEM and NOCHEM, see further details below and in 
the Appendix) 

2.1. CMIP5 Models With Interactive 
or Semi-Offline Chemistry 

[ 12 ] Eighteen of the 46 CMIP5 model versions analyzed 
here are CM1P5 models with interactive or semi-offline 
chemistry. Stratospheric and tropospheric ozone chemistry 
is interactively coupled in 9 of these 18 model versions. 
In the other nine, ozone was calculated semi-offline. We 
denote models as semi-offline if the prescribed ozone data 
set has been, unlike the models with prescribed ozone 
(see section 2.2), calculated with the underlying CMIP5 
chemistry-climate model using prescribed SSTs and SICs 
following historical emissions from Lamarque et al. [2010] 
and future emissions under the RCP scenarios as described 
by Lamarque et al. [2011]. They differ from the class of 
prescribed ozone CM1P5 models (NOCHEM), because 
their stratospheric ozone evolution responds to changes in 
GHG concentrations in the four RCPs, although it is still 
calculated offline. They are here grouped together with the 
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P 4 : Ozone is calculated with a chemistry-climate model following the four RCPs. This ozone data set is used offline in the CMIP5 model simulations. These model simulations are in the category CHEM since 
stratospheric ozone is closer to the interactive models than to the Cionni et al. [201 1] database, which follows the SRES A1B GHG scenario in the stratosphere. 

P 5 : Ozone prescribed using the data set described in Kawase et al. [2011] with tropospheric ozone following the RCPs but the impact of future climate change on stratospheric ozone not taken into account. This 
model is in category NOCHEM. 
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Total Column Ozone 




150 175 200 225 250 275 300 325 350 375 400 425 450 475 


Figure 1. Total column ozone climatologies (1986 to 2005) for (a) TOMS, (b) the N1WA combined total 
column ozone database, (c) the multimodel mean of CMIP5 models with prescribed ozone (NOCHEM), 
(d) the multimodel mean of CMIP5 models with interactive or semi-offline chemistry (CHEM), (e) the 
CCMVal-2 multimodel mean, and (f) the 1GAC/SPARC ozone database. The values in the titles are 
calculated by interpolating the models to the same grid as the TOMS observations and by subtracting 
the grid cells without observations, then calculating the area-weighted average. Individual CMIP5 models 
are shown in Figure A1 of the Appendix. 


interactive chemistry models, because the focus of the 
study is not on examining the difference between interactive 
and prescribed ozone but the sensitivity of stratospheric 
ozone evolution to GHG scenarios. The chemistry schemes 
of the interactive and semi-offline models are described in 
the Appendix (section Al.l and Table 1). 

2.2. CMIP5 Model Versions With Prescribed Ozone 

[ 13 ] In 28 of the 46 CM1P5 model versions analyzed here, 
ozone concentrations were prescribed with an offline data set 
rather than calculated interactively or semi-interactively. The 
CM1P5 models with prescribed ozone that fall into this 
NOCHEM group are driven with the four RCPs in the tropo- 
sphere. However, in contrast to the interactive and semi- 
offline models discussed in section 2.1, in the stratosphere 
the offline ozone concentrations are calculated only for one 
GHG scenario (mostly the SRES A1B). The majority of 
the models used the original or slightly modified IGAC/ 
SPARC ozone database of Cionni et al. [2011], although 
five models prescribed a different data set — the IGAC/ 
SPARC ozone data set is a zonal mean field in the 
stratosphere but three-dimensional in the troposphere (see 
Appendix section A1.2 and Table 1). 

3. Long-Term Evolution of Ozone in the 
CMIP5 Models 

[ 14 ] We start the analysis by comparing the simulated ozone 
over the past decades to observations (section 3.1) and then 
discuss projections during the 21st century in section 3.2. 


3.1. Evaluation of Ozone in the CMIP5 
Historical Simulations 

3.1.1. Evaluation of Stratospheric Ozone 

[ 15 ] Figure 1 compares the climatological mean total 
column ozone from 1986 to 2005 from the CHEM and 
NOCHEM multimodel mean against the National Institute 
of Water and Atmospheric Research (NIWA) assimilated 
data set \Bodeker et al., 2005] and observations from the 
Total Ozone Mapping Spectrometer (TOMS) [Stolarski 
and Frith, 2006]. We divide the models into these two 
classes since the ozone projections shown in section 3.2 
are done only with the CHEM models. The NIWA database 
includes measurements from four TOMS instruments, 
Global Ozone Monitoring Experiment (GOME), four solar 
backscatter ultraviolet [SBUV(/2)j instruments, and Ozone 
Monitoring Instrument (OMI). In addition, this figure 
includes a comparison to the CCMVal-2 multimodel mean 
[Eyring et al., 2010b] and the IG AC/SPARC ozone database 
[Cionni et al., 2011], The individual CMIP5 models are 
shown in Figure Al of the Appendix. Total column ozone 
has been calculated from the 3-D CMIP5 ozone mixing ratio 
fields (tro3) and the surface pressure (ps). 

[ 16 ] Both the multimodel means of CHEM and NOCHEM 
represent the observed key features in the annual cycle well, 
with highest total column ozone values in northern spring, 
a small seasonal cycle and low ozone values in the tropics, a 
relative maximum in late winter/early spring in the southern 
midlatitudes, and the ozone hole over Antarctica in spring. In 
the NOCHEM group, not all models used the IGAC/SPARC 
database to prescribe ozone (see Appendix section A 1.2), 
which explains the differences between these two panels. 


5034 



EYRING ET AL.: OZONE AND ASSOCIATED CLIMATE IMPACTS 


The amplitude and phase of the annual cycle in individual lat- 
itude bands is also well represented, but there are some quan- 
titative differences, in particular over Antarctica in spring. 
Deviations from the NIWA observations over polar regions 
in spring are larger in the NOCHEM multimodel mean 
and the IGAC/SPARC ozone database than in the CHEM 
and CCMVal-2 multimodel means. The magnitude of the 
ozone hole is slightly underestimated in the prescribed ozone 
data set, as was reported by Cionni et al. [2011]. The IGAC/ 
SPARC database is based on observations and a regression fit 
that includes decadal trends, solar cycle, and quasi-biennial 
oscillation tenns to generate continuous fields from Randel 
and Wu [2007], and over Antarctica relies on ozonesonde 
measurements from the Syowa station located at 69° S. This 
station is at the vortex edge and thus has smaller October 
mean ozone depletion than the station at the South Pole 
[Hassler et al., 2011; Solomon et al., 2005], In the CHEM 
and CCMVal-2 multimodel mean, the ozone hole over 
Antarctica in spring is well represented, with only small biases 
from the NIWA assimilated database. Larger deviations 
exist for the individual models with interactive chemistry or 
semi-offline ozone, also for other regions (see Figure Al in 
the Appendix). 

[ 17 ] F igure 2 compares the time series of total column ozone 
from 1960 to 2005 to five different observational data sets. The 
observations include ground-based measurements (updated 
from Fioletov et al. [2002]), NASA TOMS/OMI/SBUV(/2) 
merged satellite data [Stolarski and Frith, 2006], the NIWA 
combined total column ozone database [Bodeker et al., 
2005], solar backscatter ultraviolet (SBUV, SBUV/2) 
retrievals (updated from Miller et al. [2002]), and GOME/ 
Scanning Imaging Absorption Spectrometer for Atmo- 
spheric CHartographY (SCIAMACHY)/GOME-2 [Loyola 
and Coldewey-Egbers, 2012; Loyola et al., 2009]. The 
linear trend over the period with highest ozone depletion 
(1980-2000) was calculated for models and observations 
and is given in Table 2 (note that GOME/SCIA/GOME-2 
is not included in the trend calculation because it only started 
in 1995). The comparison is shown for the annual mean in 
the tropics and southern and northern midlatitudes, for March 
in the Arctic, and October in the Antarctic. To get insight 
into the vertical distribution, we additionally show compar- 
isons of climatological mean vertical profiles at different 
latitudes (October 80°S, March 0°, March 80°N) and latitu- 
dinal profiles at 50hPa in March and October in Figure 3, 
similar to Figure Al of Evring et al. [2006], The climatol- 
ogy covers the same period used to calculate the linear 
trends (i.e., 1980-2000). 

[is] In general, the CHEM and NOCHEM multimodel 
means agree well with observations, both in terms of 
absolute values and trends (see Figure 2 and Table 2). In 
agreement with observations, smallest trends are simulated 
over the tropics and largest over Antarctica in spring. 

[ 19 ] 1. In the tropics (Figure 2b), the trend and absolute 
values of the CHEM multimodel mean agree well with 
the observations, whereas the NOCHEM multimodel mean 
and the IGAC/SPARC ozone database slightly underesti- 
mate observed absolute values. All three simulate a small 
1980-2000 linear trend, comparable to observations (see 
Table 2 for trends and associated uncertainties). The three 
data sets broadly follow the observed small modulation of 
total column ozone introduced by the 1 1 year solar cycle. 


In the two models with semi-offline ozone (CCSM4, 
IPSL-CM5-LR), this modulation by the solar cycle is miss- 
ing, since the time series are prescribed from decadal means 
(see section 2.1.2). Significant differences exist among the 
individual models, with three of the CHEM models 
(CNRM-CM5, MIROC-ESM, and GFDL-CM3) and one 
NOCHEM model (MIROC4h) deviating more than around 
20 Dobson units (DU) from the observations. As noted 
above, total column ozone has been calculated from the 
tro3 fields and could slightly differ from the model’s native 
total column ozone. In the MIROC-ESM, for example, total 
column ozone is around 5 DU lower if calculated from the 
native grid, which could be due to the steep vertical gradi- 
ent of ozone mixing ratio in the tropical lower stratosphere 
that is not captured by the coarse sampling of the tro3 
fields. The overestimation in the MIROC4h model results 
from deviations in the prescribed State University of New 
York at Albany (SUNYA) data set [Wang et al., 1995] 
which was scaled by the observed stratospheric ozone trend 
from Randel and Wu [2007] for the past. The biases in the 
CNRM-CM5 could be related to the coarse vertical resolu- 
tion, the low top of the model, and/or the simplified chem- 
istry scheme (see Appendix). The CNRM-CM5 model has 
similar large biases in ozone in other regions and altitudes. 
We therefore exclude it from the multimodel mean for 
ozone projections and do not further explain deviations. 
All models except CNRM-CM5 simulate a small negative 
trend (—0.1 to — 0.36DU/yr, see Table 2) comparable to 
that in the observations (—0.1 ±0.11 to —0.19 ± 0.12 DU/ 
yr). The biases in the individual models can be attributed 
to different altitude regions in Figure 3. In the tropical 
lower stratosphere between 100 and 30hPa, 9 out of 11 
models agree well with Halogen Occultation Experiment 
(HALOE) observations and lie within the interannual stan- 
dard deviation (lcr) of the HALOE mean (Figures 3b, 3d, 
and 3e). The GFDL-CM3 model’s high bias (10 DU) seen 
in Figure 2 results from an overestimation of the ozone 
maximum at around lOhPa. At altitudes above 7hPa, 
where the photochemical lifetime of ozone is between 5 
and 10 days and the primary odd-oxygen catalytic loss 
cycle is via NO x (e.g., Lary [1997]), the three high-top 
models (CESMl(WACCM), GFDL-CM3, and M1ROC- 
ESM-CHEM) agree well with HALOE. 

[ 20 ] 2. In the southern midlatitudes (Figure 2d), the biases 
in total column ozone for the individual models discussed 
above remain the same, with the same three models 
(CNRM-CM5, GFDL-CM3, MIROC-ESM-CHEM) being 
biased compared to observations. The excess ozone at south- 
ern midlatitudes in GFDL-CM3 is likely a result of the 
excessive strength of the Brewer-Dobson circulation from 
the tropics (where ozone is already too high) to midlatitudes. 
The CHEM multimodel mean starts at too high total column 
ozone values in 1960, similar to the CCMVal-2 multimodel 
mean. Overall, in using the Randel and Wu [2007] observa- 
tions for this period, the NOCHEM multimodel mean and 
the IGAC/SPARC ozone database are in agreement with 
the observations shown here and agree better with the obser- 
vations than the CCMVal-2 and CHEM multimodel means. 

[ 21 ] 3. In the northern midlatitudes (Figure 2c), the 
MIROC-ESM-CHEM is still biased high, but the CNRM- 
CM5 and the GFDL-CM3 are closer to observations. In this 
region, CCSM4 simulates a positive trend (0.07 ± 0.01 DU/ 


5035 


EYRING ET AL.: OZONE AND ASSOCIATED CLIMATE IMPACTS 


320 


310 


300 

D 

O 

290 


280 


270 



260 

1960 1970 1980 1990 2000 


330 


Total Ozone Column ANN 90N-90S 


Year 


380 

360 

Q 3^0 

320 

300 



Total Ozone Column ANN 60N-35N 


280 

I960 1970 1980 1990 2000 

Year 


550 
500 
Q “50 
400 
350 



600 


Total Ozone Column MAR 90N-60N 


300 

I960 1970 1980 

Year 


1990 2000 


Total Ozone Column ANN 25N-25S 



240 - - "" 

220 - 

200 


1960 1970 1980 1990 2000 

Year 


380 
360 
Q 340 
320 
300 



280 

1960 1970 1980 1990 2000 

Year 


400 


Total Ozone Column ANN 35S-60S 


450 


400 


350 

Z) 

o 

300 


250 


200 



150 

1960 1970 1980 1990 2000 

Year 


500 


Total Ozone Column OCT 60S-90S 


IGAC/SPARC CMIP5Chem 

' CMIP5noChem MIROC4h MIROC-ESM-CHEM 

GISS-E2-R.p3 GISS-E2-R.p2 

GISS-E2-H.p3 GISS-E2-H.p2 GFDL-CM3 

CNRM-CM5 CESM1-FASTCHEM CESM1-WACCM 

- CCSM4 * DLR_GOME-SCIA-GOME2 V NIWA 


□ NASATOMS-SBUV-OMI o NOAA_SBUV + Ground-based 


Figure 2. Time series of area-weighted total column ozone from 1960 to 2005 for the annual mean 
(a) global mean (90°S-90°N), (b) in the tropics (25°S-25°N), (c) in the northern midlatitudes (35°N-60°N), 
(d) in the southern midlatitudes (35°-60°S), and for the March and October mean in the (e) Arctic 
(60°N-90°N), and (f) Antarctic (60°S-90°S), respectively. The figure shows the multimodel mean and 
individual CM1P5 models with interactive or semi-offline chemistry (CHEM, red solid and colored lines) 
and standard deviation (blue shaded area) compared to the IGAC/SPARC ozone database (black dotted 
line). In addition, the multimodel mean of the CMIP5 models that prescribe ozone (NOCHEM, green 
solid line), the CCMVal-2 multimodel mean (yellow solid line), and observations from five different 
sources (see text) are shown. 
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Table 2. Linear Trends and Errors in Total Column Ozone (DU/year) d 



90°S-90°N 

AM 

25°N-25°S Annual 
Mean 

35°N-60°N Annual 
Mean 

60°S-35°S Annual 
Mean 

60°N-90°N March 
mean 

90°S-60°S Oct 
mean 

Ground-based 

-0.56 ±0.11 

-0.08 ±0.10 

-0.84 ±0.25 

-0.90 ±0.17 

-3.18 ±0.95 

-5.25 ±0.77 

NOAA-SBUV 

-0.74 ±0.12 

0.16 ±0.13 

1.12±0.21 

-1.21 ±0.18 

-3.30 ±0.96 

-3.99 ±0.80 

NASA TOMS-SBUV- 
OMI 

-0.67 ±0.13 

-0.19 ±0.12 

-0.89 ±0.23 

-1.06 ±0.17 

-2.95 ±0.91 

-3.70 ±0.82 

NIWA 

0.61 ±0.12 

0.10 ±0.11 

-0.88 ±0.23 

-0.87 ±0.16 

-3.18 ±0.92 

-3.80 ±0.78 

CCSM4 

— 0. 1 7 ± 0.02 

-0.08 ±0.00 

0.07 ± 0.01 

-0.35 ±0.03 

-0.60 ±0.04 

-2.63 ±0.26 

CESM 1 (WACCM) 

-0.52 ±0.05 

-0.13 ±0.06 

-0.31 ±0.16 

-0.94 ±0.12 

— 1.85 ± 0.86 

— 5.60± 1.25 

CESMl(FASTCHEM) 

— 0.56±0.16 

-0.36 ±0.14 

-0.34 ±0.34 

-0.78 ±0.25 

— 0.74± 1.17 

—4.12 ± 0.78 

CNRM-CM5 

-0.00 ±0.03 

0.12 ±0.03 

0.16±0.09 

— 0.19±0.06 

-0.03 ±0.38 

— 2.30± 0.37 

GFDL-CM3 

-0.79 ±0.07 

-0.31 ±0.07 

—0.52 ±0.15 

— 1.46 ± 0. 14 

— 1.78 ± 0.94 

-7.33 ±0.98 

GISS-E2-H.p2 

— 1 .05 ±0.17 

-0.33 ±0.10 

-1.21 ±0.27 

-1.86 ±0.22 

-3.28 ± 1.04 

— 5.34± 0.72 

GISS-E2-H.p3 

— 1 .06 ±0.19 

-0.36 ±0.09 

-1.06 ±0.34 

-2.04 ±0.30 

-2.96 ±0.99 

— 5.54± 0.86 

GISS-E2-R.p2 

— 1.11 ±0.19 

-0.33 ±0.08 

-1.28 ±0.32 

-2.07 ±0.30 

-2.84 ± 1.08 

-6.08 ±0.91 

GISS-E2-R.p3 

— 1.18 ± 0.19 

-0.32 ±0.08 

-1.33 ±0.34 

-2.28 ±0.27 

—3.27 ± 1.32 

-7.07 ±0.91 

IPSL-CM5-ALR 

-0.43 ± 0.03 

-0.10±0.01 

— 0. 18 ± 0.02 

-0.96 ±0.06 

— 1.08 ± 0.07 

-5.56 ±0.32 

MIROC-ESM-CHEM 

-0.55 ±0.07 

-0.17 ±0.07 

—0.37 ± 0.1 1 

— 1.1 7 ±0.15 

— 0.80±0.58 

— 5.24± 0.70 

MIROC4h 

-0.82 ±0.04 

-0.05 ±0.02 

-1.26 ±0.04 

-1.06 ±0.04 

—3.85 ± 0.14 

-6.51 ±0.21 

NOCHEM 

-0.63 ±0.05 

-0.22 ±0.04 

-0.84 ±0.08 

-0.83 ±0.04 

-2.00 ±0.08 

-4.63 ±0.21 

CHEM 

-0.67 ±0.09 

-0.21 ±0.04 

— 0.58±0.13 

— 1.28 ± 0.13 

— 1.75 ± 0.41 

— 5.16 ± 0.38 

CCMVal-2 

-0.61 ±0.03 

-0.32 ±0.03 

-0.56 ±0.04 

-0.97 ±0.05 

— 1.45 ± 0.21 

-4.08 ±0.26 

IGAC/SPARC 

-0.64 ±0.06 

-0.33 ±0.05 

-0.79 ±0.09 

-0.83 ±0.04 

— 1 .57 ± 0.08 

-4.53 ±0.22 


a Over the 1980-2000 period for observations, the individual CMIP5 models, the CHEM, NOCHEM, and CCMVal-2 multimodel means, and the IGAC/ 
SPARC database for the different regions shown in Figure 2. 


yr) over the 1980 to 2000 period in contrast to the observations 
which show a negative trend (—1.12 ± 0.21 to —0.84 ± 0.25 
DU/yr) and is also biased low in tenns of absolute values by 
around 20 DU at the beginning of the time series in the 
1960s. This is partially related to the fact that the model top 
is located at 2.2 hPa. There is also the factor that the fields in 
CCSM4 are linearly interpolated in time from decadal aver- 
ages from the NCAR-CAM3.5 results, smoothing out the re- 
sults (see section 2. 1.2.1). 

[ 22 ] 4. Over Antarctica in October (Figure 2f), CCSM4 
and CNRM-CM5 substantially underestimate ozone deple- 
tion (—2.63 ±0.26 and —2.3 ± 0.37 DU/yr, respectively, 
compared to observed -5.25 ±0.77 to -3.70 ± 0.82 DU/ 
yr). This underestimation will lead to subsequent biases in 
zonal mean wind changes [Son et al., 2008] which are 
further discussed in section 4. CESMl(FASTCHEM) is 
biased slightly low, in particular at the beginning of the time 
series in the 1960s. This is presumably because it uses a simple 
temperature threshold and decay to represent ozone hole loss 
processes [Cariolle et al., 1990] so the ozone loss starts in 
the winter before the sun rises. All other models and the 
CHEM, NOCHEM, and CCMVal-2 multimodel means agree 
well with observations, both in terms of trends and 
absolute values, although the CMIP5 CHEM multimodel 
mean represents the ozone depletion period slightly better 
(see also Figure 1). Ozone recovery starts too early in the 
CCSM4 and 1PSL in around 1995. These are the two models 
where interannual variability in ozone is missing, since the 
data sets stem from decadal means. This could lead to biases 
in temperature trends as is further discussed below. 

[ 23 ] 5 . hi the Arctic in March (Figure 2e), again CCSM4 and 
CNRM-CM5 underestimate ozone depletion (—0.60 ± 0.04 
and —0.03 ± 0.38, respectively, compared to observed 
—3.30 ± 0.96 to —2.95 ± 0.91 DU/yr). All other models also 
simulate smaller trends than observed, leading to an 
underestimation of the trend by the CHEM multimodel mean 
(—1.75 ± 0.41 DU/yr) and the NOCHEM multimodel mean 


(—2.00 ± 0.08 DU/yr), and also a similar underestimation 
by the IGAC/SPARC ozone database (-1.57 ± 0.08 DU/ 
yr). The absolute values are mostly of the order of those ob- 
served, except for CESMl(WACCM) and GFDL-CM3, 
which are biased low. 

[ 24 ] 6. The global annual mean (Figure 2a) combines the 
different regions discussed above. In terms of absolute 
values, the CHEM and CCMVal-2 multimodel means agree 
well with observations. The NOCHEM multimodel mean 
and the IGAC/SPARC ozone database are biased slightly 
low, by around 10 DU. The 1980-2000 trends of all four time 
series lie between -0.61 ±0.03 and -0.67 ± 0.09 DU/yr, 
within the observed range (—0.56 ±0.11 to —0.74 ±0.12 
DU/yr). The models also consistently simulate a negative 
trend (—0.17 ±0.02 to —1.18 ± 0.19 DU/yr) except for 
CNRM-CM5 that has a zero trend (±0.03). 

3.1.2. Evaluation of Tropospheric Ozone 

[ 25 ] As mentioned, ozone in both the CHEM and 
NOCHEM group follows the RCP scenarios in the future. 
The tropospheric ozone portion of the IGAC/SPARC 
ozone database, which was used in the majority of the models 
to prescribe ozone (see Table 1), is based on calculations of 
the historical period from CAM3.5 [Lamarque et al., 2010] 
and the NASA GISS model for Physical Understanding 
of Composition-Climate Interactions and Impacts (GISS- 
PUCCINI, Shindell et al. [2006]); see further description 
by Cionni et al. [2011], The three-dimensional (latitude, 
longitude, altitude) future tropospheric ozone time series 
continues the historical CAM3.5 simulation (not GISS- 
PUCCINI) until 2099, following the four different RCPs. 
To be consistent with the ACCM1P calculations of Young 
et al. [2013b], we used the WMO tropopause definition to 
identify the tropopause. Here, the boundary between the 
troposphere and the stratosphere is defined as the lowest level 
at which the lapse rate decreases to 2 K/km or less, provided 
that the average lapse rate between this level and all higher 
levels within 2 km does not exceed 2 K/km. Figure 4 shows 
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Figure 3. Climatological zonal mean (1980-2000) ozone mixing ratios from the individual CMIP5 
models, the CHEM and NOCHEM multimodel means, and the IGAC/SPARC database compared 
to HALOE observations in pmol/mol. Vertical profiles at (a) 80°N in March, (b) 0° in March, and 
(c) 80° S in October. Latitudinal profiles at 50hPa in (d) March and (e) October. The gray area shows 
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a comparison of simulated tropospheric column ozone values 
averaged over the period 2000-2005 to values derived from 
OMI and Microwave Limb Sounder (MLS) instruments on 
board the Aura satellite averaged over 2000-2005 ( Ziemke 
et al. [2011]; source: http://acdb-ext.gsfc.nasa.gov/Data ser- 
vices/cloud slice/new data.html). The individual models are 
shown in Figure A2 of the Appendix. 

[ 26 ] Figure 4 shows that the geographical distribution of 
tropospheric column ozone is very similar in the CHEM 
and NOCHEM multimodel means. They agree in important 
aspects with the observations, for example, have lower values 
in the tropics and lower values in the SH compared to the 
Northern Hemisphere (NH). Near-global (60°S-60°N) mean 
absolute values nearly agree, but this is to some degree due 
to the compensation of a high model bias in NH midlatitudes 
and a low model bias in the tropics. The models overall 
show a stronger hemispheric gradient in tropospheric ozone 
than is seen in the observations. As with comparisons against 


stratospheric ozone data, substantial differences exist among 
the individual models and compared to observations (see 
Figure A2 in Appendix). Most notably, the 1PSL models use 
zonal mean ozone in the troposphere and therefore miss the 
observed longitudinal structure, such as ozone maxima over 
regions with high ozone precursor emissions and the very 
low ozone over the tropical Pacific Ocean which has been 
previously detected (e.g., Singh et al. [1996]). The spatial 
distribution of ozone could be important for estimates of 
regional radiative forcing (RF) [Freckleton et al., 1998; 
Shindell and Faluvegi, 2009] and OH levels and hence life- 
times for many very short lived species (VSLSs). In addition, 
the differences in absolute ozone values among the CMIP5 
models could lead to differences in global mean tropospheric 
ozone RF, discussed further by Young et al. [2013b]. 

[ 27 ] Figure 5 compares a 20 year (1986-2005) climatol- 
ogy of the CHEM models and the IGAC/SPARC ozone 
database against ozonesonde measurements, which broadly 


5038 


EYRING ET AL.: OZONE AND ASSOCIATED CLIMATE IMPACTS 



CMIP5noChem 28.2 DU 

90N 


CMIP5Chem 30.5 DU 


60N 

SON 

0 

SOS 

60S 


1 80 1 50WI 20W90W 60W SOW 0 30E 60S 90S 1 20E1 50E 1 80 


180150W20W90W60W30W 0 30E 60E 90E 120E150E 180 


90N 

SON 

SON 

0 

30S 

60S 



Figure 4. Climatological mean annual mean tropospheric column ozone averaged between 2000 and 
2005 from the CMIP5 NOCHEM and the CHEM multimodel mean compared to the 2005-2009 MLS/ 
OMI observations [from Ziemke et al., 2011], The numbers in the title show the near-global area- 
weighted mean (60°S-60°N) as in the observations. Individual CMIP5 models are shown in Figure A2 
of the Appendix. 


cover the same time period and are described by Logan 
[1999], The models were sampled at the ozonesonde sites, 
and the observational and model data were aggregated into 
four latitude bands; data from the lower (750 hPa), mid 
(500 hPa), and upper (250 hPa) troposphere are shown. A 
similar summary figure was presented by Stevenson et al. 
[2006] and Young et al. [2013b]. The multimodel mean 
captures the seasonal cycle in observed ozone concentrations 
well, being significantly correlated (r> 0.58; 5% level) with 
the ozonesondes in all locations except the NH tropical 
upper troposphere where the correlation is only 0.45. Young 
et al. [2013b] showed that the multimodel mean from the 15 
state-of-the-art chemistry models participating in ACCMIP 
also could not capture the seasonal cycle here, perhaps 
indicating problems with the modeling of the relevant 
processes for this region (e.g., convective transport, lightning 
NO x , natural emissions, and stratosphere-troposphere connec- 
tions), although we note that the interannual variation is 
comparatively large in this region. Compared to the climato- 
logical mean for the observations, there are negative biases 
for the multimodel mean in the SH and mostly positive biases 
in the NH (except the tropical upper troposphere again), which 
is also in agreement with Young et al. [2013b]. However, we 
note that the multimodel mean is within the interannual vari- 
ability of the majority of the observations, except for winter 
and spring in the NH extratropics. The spread in 
performance of the individual CMIP5 models is generally 
larger than in similar multimodel comparisons [Stevenson 
et al., 2006; Young et ah, 2013b], with CESM 1 (FASTCHEM) 
and CNRM-CM5 notably underestimating the observations 
through much of the tropics and the GISS models notably 
overestimating the observations in the NH extratropics. For 
CESM1 (FASTCHEM), this is probably because lightning 
NO x was inadvertently missing in these simulations (see 


section 2.1.1. 1) and differences in lightning NO x or isoprene 
emission could also be possible reasons for the deviations in 
the GISS model. Ozone concentrations from the IGAC/ 
SPARC ozone database are lower than the multimodel mean 
and fall outside the range of concentrations described by 
interannual variability in the SH tropical upper troposphere. 

3.2. Future Evolution of Ozone 

3.2.1. Stratospheric Ozone Projections 

[ 28 ] The future evolution of stratospheric ozone is 
assessed in the concentration driven historical and RCP 
simulations of the models with interactive or semi-offline 
chemistry. However, we exclude the CNRM-CM5 model 
in this group, since the simplified chemistry scheme leads 
to poor performance of ozone compared to observations 
(see section 3.1). We also exclude the IPSL-CM5A model, 
since it uses a simplified scaling to determine ozone for the 
various RCPs (see section Al.2.2). In contrast to the 
historical simulations, future changes in tropospheric column 
ozone contribute a significant amount to changes in total 
column ozone under the RCPs (see section 3.2.2). Conse- 
quently, we show stratospheric column instead of total column 
ozone used in the previous section. In the evaluation section, 
total column ozone was chosen since long-term observations 
exist for this data product. 

[ 29 ] The evolution of stratospheric ozone over the histori- 
cal time period (1960 to 2005) and the sensitivity of ozone to 
future GHGs (2006-2100) for the four different RCPs are 
first analyzed for the six regions discussed above (Figure 6). 
We note again that all simulations follow the same ODS sce- 
nario in the future, namely the Al scenario of WMO [2003], 
The CHEM models that are included in the different RCPs 
for future analysis are listed in the caption of Figure 6, and 
resulting historical and future changes among the RCPs for 
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Figure 5. Comparison of the annual cycle of ozone, between observations (black circles) and the CMIP5 
multimodel mean (solid red line) and CMIP5 multimodel median (dashed red line) averaged over the period 
1986-2005. Model and observational data were grouped into four latitude bands (90°S to 30°S, 30°S to 0°, 
0° to 30°N, and 30°N to 90°N) and sampled at three altitudes (700 hPa, 500 hPa, and 250 hPa), with the 
models sampled at the ozonesonde locations before averaging together. The individual CHEM models are 
represented by the colored solid and dashed lines. Error bars on the observations indicate the average 
interannual standard deviation for each group of stations. The correlation (r) and mean normalized bias error 
(mnbe) for the CMIP5 multimodel mean versus the observations are also indicated in each panel. This figure 
is similar to F igure 4 of Young et al. [20 1 3b], who discuss the ACCMIP instead of the CM1P5 CHEM models. 


the corresponding multimodel global annual mean are 
summarized in Table 3. 

[ 30 ] To calculate multimodel trends and their confidence 
and prediction uncertainties, we use the same time series 
additive model (TSAM) method [Scinocca et al., 2010] that 
was applied in previous CCMVal studies [Eyring et al., 
2010b; SPARC-CCMVal, 2010], We also use the same 
terminology, where the term “trend” refers to a smooth 
trajectory passing through the model time series. This repre- 
sents the signal resulting from forced changes, leaving noise 
as a residual resulting from internal, unforced climate vari- 
ability. Since we are interested in the question whether the 
long-term evolution differs between the RCPs, we show 
the point-wise 95% confidence interval. This interval has a 
95% chance of overlapping the true trend, representing the 
local uncertainty in the trend at each year. The TSAM 
method was also previously applied in a study on the sensi- 
tivity of stratospheric ozone to GHG scenarios but applying 


a smaller set of models [Eyring et al., 2010a]. We reassess 
their conclusions with the help of the CM1P5 models with 
interactive and semi-offline chemistry in comparison to 
the IGAC/SPARC ozone database, which was used in the 
majority of the CM1P5 models with prescribed ozone. It 
should be noted that the model used to calculate semi- 
offline ozone in CCSM4 (CAM3.5) was also included 
in Eyring et al. [2010a]. Similar to section 3.1, we 
additionally show vertical profiles (here: 2090s-2000s and 
2090s-1960s) to probe the vertical distribution of future 
stratospheric ozone changes. 

[ 31 ] There are no significant changes prior to 1960 
(not shown). From 1960 onwards, stratospheric ozone 
decreases due to increase in ODSs until around 2000, except 
in the tropics where changes are generally small. This 
decrease is followed by a steady increase with a rate being 
dependent on the RCP and the region. Figure 7 shows the 
differences in vertically resolved ozone between the 2090s 
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Figure 6. The 1980 baseline-adjusted stratospheric column ozone time series from 1960 to 2100 for the 
CM1P5 CHEM multimodel mean (colored lines) and the IGAC/SPARC ozone database (black line) for dif- 
ferent regions. All time series by construction go through 0 in 1980. The RCP 2.6 (CCSM4, GFDL-CM3, 
GISS-E2-R.p2, GISS-E2-R.p3, M1ROC-ESM-CHEM), RCP 4.5 (CCSM4, CESM1-WACCM, GFDL- 
CM3, GISS-E2-R.p2, GISS-E2-R.p3, MIROC-ESM-CHEM), 6.0 (CCSM4, GFDL-CM3, GISS-E2-R.p2, 
GISS-E2-R.p3, MIROC-ESM-CHEM), and RCP 8.5 (CCSM4, CESM1-WACCM, GFDL-CM3, GISS- 
E2-R.p2, GISS-E2-R.p3, MIROC-ESM-CHEM) are shown in blue, light blue, orange, and red, respectively, 
with the models included in the multimodel mean being shown in brackets. The corresponding color-coded 
stippled areas show the 95% confidence interval of the CHEM multimodel mean simulations. 


and the pre-ozone depleted 1960s (left two columns) and 
between the 2090s and the 2000s (right two columns). 

[ 32 ] 1. In all RCP simulations, a decrease in the tropical 
lower stratospheric ozone is simulated in the altitude range 
between around 100 and 30hPa compared to both the 
2000s (Figure 7, left) and the 1960s (Figure 7, right). This 
decrease is smaller in the two lower (RCP 2.6 and 4.5) 
GHG scenarios compared to the higher ones (RCP 6.0 and 
8.5). The decrease in the lower tropical stratosphere is a 


result of a projected strengthening of the Brewer-Dobson 
circulation [ Butchart et al., 2010; Butchart et al., 2011; 
SPARC-CCMVal, 2010]. In contrast, in the upper strato- 
sphere the sensitivity to GHG scenarios is opposite in sign, 
with smaller GHG-induced cooling leading to smaller 
changes in ozone from changes in ozone destruction rates. 
Compared to pre-ozone depletion (1960s), this increase 
quantifies “super-recovery,” which is a C0 2 -induced cooling 
effect in the upper and middle stratosphere that slows down 
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Table 3. Stratospheric Column Ozone (DU) a 


Year 

RCP 2.6 

RCP 4.5 

RCP 6.0 

RCP 8.5 

1850 

293 ±9 

293 ±9 

293 ±9 

293 ±9 

1980 

291 ±8 

291 ±8 

291 ±8 

291 ±8 

2000 

276 ±8 

276 ±8 

276 ±8 

276 ±8 

2005 

276 ±9 

276 ±9 

276 ±9 

276 ±9 

2010 

278 ±9 

275 ± 1 1 

277 ±9 

274 ± 12 

2020 

280 ± 10 

276 ± 12 

279 ± 1 1 

278 ± 11 

2030 

284 ±9 

280 ± 10 

283 ±9 

281 ±10 

2040 

285 ±9 

283 ± 1 1 

286 ± 10 

285 ± 10 

2050 

288 ±9 

286 ± 12 

289 ± 10 

288 ± 11 

2060 

289 ±9 

288 ± 11 

290 ±9 

291 ±11 

2070 

289 ±9 

287 ±9 

291 ±9 

292 ±8 

2080 

288 ±9 

287± 10 

292 ±9 

291 ±11 

2090 

289 ±9 

288± 10 

292 ± 10 

294 ± 10 

2100 

291 ±9 

289 ± 11 

293 ± 10 

296 ±11 


a For the various RCPs as calculated from the multimodel mean of the 
CMIP5 models with interactive or semi-offline stratospheric chemistry 
(see caption of Figure 6). The multimodel mean is given along with uncer- 
tainties (±1 SD). 

ozone destruction rates, leading to an increase in ozone. The 
resulting evolution of tropical stratospheric ozone column 
that can be seen in Figure 6 is a combination of decreases 
in lower stratospheric ozone and increases in upper strato- 
spheric ozone that result in small changes in the total column 
over the 21st century. The decrease that was projected by the 
CCMVal-2 models due to an enhanced Brewer-Dobson cir- 
culation under the SRES A1B scenario [SPARC-CCMVal, 
2010] is again seen in the higher RCP 6.0 and 8.5 scenarios, 
in which CO 2 concentrations continue to increase in the 
second half of the 21st century (see for example Figure 1 
of Cionni et al. [2011]). This decrease is not simulated in 
the RCP 4.5 and RCP 2.6 scenarios, where CO 2 concentra- 
tions level or slightly decrease in the second half of the 
21st century, respectively. Plowever, notably also in the 
RCP 2.6 and 4.5 scenarios, tropical ozone remains below 
1960 levels at the end of the 21st century. These findings 
agree well with current general understanding on processes 
in the upper and lower stratosphere [Butchart et al., 2010; 
Eyring et al., 2010a; SPARC-CCMVal, 2010]. 

[ 33 ] 2. In the midlatitudes, the ozone sensitivity to GHG 
scenarios in the upper stratosphere is the same as in the tro- 
pics. The upper stratospheric behavior is again related to 
stronger C02-induced cooling in the higher GF1G scenarios, 
in line with previous findings [SPARC-CCMVal, 2010; 
WMO, 2011], In contrast to the tropics, in the lower strato- 
sphere substantial increases compared to the 2000s (Figure 7, 
left) and the 1960s (Figure 7, right) are simulated. These 
increases are higher in the scenarios with higher GHG con- 
centrations. The substantial increase of UTLS ozone as seen 
particularly in the higher RCP scenarios is important since 
ozone changes in this altitude region have the greatest impact 
on surface temperature on purely radiative grounds [Forster 
and Shine, 1997]. The increase could again be due to an 
increased Brewer-Dobson circulation in the higher GHG 
scenarios, resulting in an increased stratosphere-troposphere 
exchange and ozone flux. This increase is smaller in the 
southern than in the northern midlatitudes [Shepherd, 2008], 
Overall, the positive ozone sensitivity to GHGs in the upper 
and lower stratosphere leads to an increase in stratospheric 
column ozone that is positively correlated with the GHG 
scenarios (Figure 6c). 


[ 34 ] 3 . Over the Arctic in March, the results are very sim- 
ilar to the northern midlatitudes, with a more pronounced 
sensitivity resulting in differences by 2100 that are larger 
than in the northern midlatitudes (~40 DU among the RCPs). 
Despite this difference, the ozone return dates are mostly 
similar among the RCPs, since stratospheric ozone returns 
to pre-1980 values before the time when the scenarios 
deviate. Substantial super-recovery is simulated in the higher 
RCP scenarios. The IGAC/SPARC ozone database is overall 
very similar to RCP 6.0. 

[ 35 ] 4. Over Antarctica in September, the sensitivity of 
ozone to the GHG scenario is generally not very high, 
resulting in very similar evolution of stratospheric column 
ozone over the 21st century in the four RCPs, with ozone 
starting to increase over the coming decades and returning 
back to pre-1960 values only by the end of the 21st century, 
as has been found in previous studies [Austin et al, 2010; 
SPARC-CCMVal, 2010; WMO, 2011], Similar to the 
midlatitudes and the Arctic, total column ozone values are 
highest in the RCP 8.5 at the end of the century. 
Differences among the RCPs until 2050 are very small, 
confirming that differences in GHGs are more important in 
the second half of the 21st century when the impact of ODSs 
due to the projected decrease slowly disappears. 

3.2.2. Tropospheric Column Ozone Projections 

[36] Unlike the stratospheric halogen scenarios, the 
precursor emissions are not the same in the different RCPs 
(see, e.g., Figure 1 of Cionni et al. [201 1]). While it is clear 
that by 2100 all RCPs have substantially smaller NO x , CO, 
and VOC emission totals compared to 2000, the largest 
emissions are in RCP 8.5. Furthermore, another major 
difference with RCP 8.5 is the change in the ozone precursor 
CH 4 , which nearly doubles between 2000 and 2100, while 
remaining at similar or decreased levels in the other three 
RCPs. Some of these factors are evident from comparing 
the panels in Figure 8, which presents decadal averages for 
the tropospheric ozone column in the CHEM mean for the 
1940s and 2000s (historical), and 2040s and 2090s (RCPs). 
An extended list of decadal multimodel global mean values 
is given in Table 4, including standard deviations. This 
figure is similar to Figure 9 of Cionni et al. [2011] but uses 
the larger CHEM mean in contrast to a single model projec- 
tion (CAM3.5), although the general evolution and key 
findings of that study remain. As noted previously [Cionni 
et al., 2011; Kawase et al., 2011; Lamarque et al., 2011; 
Young et al., 2013b], the large differences in tropospheric 
column ozone between RCP 8.5 and the other RCPs are 
likely due in part to the increased mixing ratios of the ozone 
precursor methane in RCP 8.5, which act to strongly 
increase tropospheric ozone globally, although the largest 
increases are around 30°N over Europe, Africa, and Far East 
Asia, which could relate to the redistribution of precursor 
emissions. In addition, the important role of an increase 
in the input of stratospheric ozone, driven both by the 
stratospheric ozone recovery [SPARC-CCMVal, 2010; 
WMO, 2011] and the strengthening in the residual circula- 
tion with increased GHGs [Butchart et al., 2010; Butchart 
et al., 2011], has been noted [Kawase et al., 2011; Young 
et al., 2013b], This process could be driving the approxi- 
mately zonal increases in the tropospheric ozone column 
at midlatitudes, where much stratospheric ozone enters the 
troposphere [Olsen et al., 2004]. For RCP 2.6, tropospheric 
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Figure 7. Differences in vertically resolved ozone between the 2090s and 2000s (Figures 7al to 7fl) and 
2090s and 1960s (Figures 7a2 to 7f2) for the RCP simulations for (al and a2) global (90°S-90°N) annual 
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mean, and (fl and f2) Antarctic (60°S-90°S) October mean. 


column ozone in the 2090s looks very similar to the 1940s, 
reflecting the strong decreases in precursor emissions 
(including methane) compared to 2000. Smaller decreases 
in precursor emissions, approximately level methane, and, 
potentially, increased stratosphere-to-troposphere transport 
of ozone cause RCP 4.5 and 6.0 tropospheric column 
ozone to remain comparable to the 2000s by 2 100. We 
have briefly discussed tropospheric column ozone projec- 
tions here, since differences in tropospheric column ozone 
(ranging from 27.8 DU in the RCP 2.6 and 43.8 DU in the 
RCP 8.5 by 2100, see Figure 8) contribute substantially 


to differences in total column ozone projections among 
the RCPs. We therefore argue that total column ozone 
should not be used to project stratospheric ozone in the 
RCPs. 

4. Climate Impacts Associated With Stratospheric 
Ozone Changes 

[ 37 ] Previous studies have built a strong case for the 
impact of stratospheric ozone changes on climate. Here, 
we highlight the role of ozone on climate in the consistent 
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Figure 8. Decadal mean tropospheric column ozone for the 1950s and 2000s (left) and the RCPs (right); 
from top to bottom: RCP 2.6, RCP 4.5, RCP 6.0, and RCP8.5 shown for the 2050s and 2090s. The data 
are derived from the CHEM multimodel mean. For the different RCPs, the same models as in Figure 6 are 
included in the CMIP5 CF1EM multimodel mean. A similar figure showing the IGAC/SPARC database 
can be found in Cionni et al. [2011]. 


framework of CM1P5, which includes for the first time 
multiple climate models that are both coupled to an ocean 
and have interactive chemistry and contrast these with the 
results of the CCMVal-2 simulations [Son et al., 2010]. We 
highlight that variations in ozone forcing can play a large role 
in some aspects of the resulting climate projections and 
caution that ozone differences need to be accounted for in 
any interpretation of associated climate shifts. 

4.1. Impact of Ozone on Lower Stratospheric 
Temperature Trends 

[38] Stratospheric ozone depletion is a key forcing in 
explaining the observed evolution of lower stratospheric global 
mean temperature during the satellite period, being the main 
contributor to the long-term cooling trend [Gillett et al., 2011; 
Ramaswamy et al., 2006]. Furthermore, observations indicate 
that the temporal variability in global mean ozone immediately 
following volcanic eruptions contributed to the step-like 
temperature decline of global mean stratospheric temperatures 
in the aftermath of the eruptions of El Chichon in 1982 and 
Mt Pinatubo in 1991 [Eyring et al., 2006; Thompson and 
Solomon, 2009], Stratospheric ozone changes could also play 
a key role in tropical lower stratospheric temperature trends 
[Polvani and Solomon, 2012; Santer et al., 2012; Solomon 
et al., 2012; Thompson et al., 2012]. 

[ 39 ] To assess the impact of ozone in the CMIP5 models 
on the climate of the recent past, we first evaluate time 
series of temperatures in the lower stratosphere (TLS) 


for near-global annual averages (Figure 9a) and October 
through January averages over SH polar latitudes 
(Figure 9b). We compare observational estimates based on 
Microwave Sounding Units (MSU) of the Remote Sensing 
Systems (RSS) version 3.3 [Means et al., 2011] with 
synthetic MSU temperatures derived from climate model 
output with the TLS vertical weighting function applied to 


Table 4. Tropospheric Column Ozone (DU) a 


Year 

RCP 2.6 

RCP 4.5 

RCP 6.0 

RCP 8.5 

1850 

21.3 ± 3.6 

21.3 ±3.6 

21.3 ± 3.6 

21.3 ± 3.6 

1980 

29.5 ±3.7 

29.5 ±3.7 

29.5 ±3.7 

29.5 ±3.7 

2000 

31.1 ±3.3 

31.1 ±3.3 

31.1 ±3.3 

31.1 ±3.3 

2005 

31.5 ± 3.5 

31.5 ±3.5 

31.5 ± 3.5 

31.5 ± 3.5 

2010 

32.2 ±2.9 

3 1 .7 ± 2.5 

31.9 ± 2.7 

31.9 ± 2.7 

2020 

32.1 ±2.6 

32.0 ±2.8 

32.1 ±2.8 

33.2 ±3.1 

2030 

31.7 ±2.9 

32.6 ±2.8 

32.5 ±3.2 

34.6 ±3.6 

2040 

31.6 ± 3.0 

32.7 ±2.9 

33.2 ±3.7 

35.6 ±3.7 

2050 

31.1 ±2.9 

32.8 ±3.2 

33.5 ±3.8 

36.8 ±3.9 

2060 

30.4 ±3.1 

32.4 ±3.3 

33.7 ± 3.9 

38.2 ±4.5 

2070 

29.5 ±2.9 

3 1.6 ± 3. 1 

33.4 ±4.2 

39.2 ±4.4 

2080 

28.6 ±2.9 

3 1 .0 ± 3.3 

33.1 ±4.2 

40.0 ±4.8 

2090 

28.3 ±3.0 

30.7 ±3.4 

32.6 ±4.0 

40.6 ±4.7 

2100 

28.0 ±2.8 

30.6 ±3.2 

32.2 ±4.2 

41.3 ±4.9 


“For the various RCPs as calculated from the multimodel mean of the 
CMIP5 models with interactive or semi-offline tropospheric chemistry 
(see caption of Figure 6). The multimodel mean is given along with uncer- 
tainties (±1 SD). 
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Figure 9. Time series from 1979-2005 of (a) near-global 
(82.5°S to 82.5°N) annual mean and (b) Southern Hemisphere 
polar (82.5°S to 60°S) October-January mean lower strato- 
spheric temperatures (TLS) for the CMIP5 multimodel mean 
compared to RSS observations (black dots). The individual 
models of the CHEM group are shown in blue colored lines, 
while those from the NOCHEM group are shown in orange. 
The time series in Figure 9b are smoothed using a 3 year 
running average. The multimodel mean time series are not 
shown due to the lack of stratospheric aerosol forcing in some 
of the model simulations. 


the near-global (82.5°S-82.5°N) domain of the MSU data. 

[ 40 ] The CM1P5 models, when also including stratospheric 
aerosol forcing, represent well the observed step-like reduction 
of global mean TLS (Figure 9a). The CM1P5 multimodel 
mean trend over the period (—0.28 ± 0.08 K/decade) underes- 
timates the long-tenn cooling trend from the RSS satellite 
observations (—0.35 ± 0.07 K/decade), consistent with Santer 
et al. [2012], On average, the linear trend is slightly smaller 
in the CHEM group (—0.26 ±0.11 K/decade) than in the 
NOCHEM group (—0.30 ± 0.05 K/decade). We note, how- 
ever, that several models exhibit trend values similar to or 
larger than the observed RSS TLS estimate. The wide range 
of simulated global TLS historical trends is consistent with a 
large spread of simulated ozone trends (Table 2). This is 
depicted in Figures 10a and 10b, which examines the 
multimodel mean near-global (82.5°S-82.5°N) lower strato- 
spheric ozone (50hPa) and temperature (lOOhPa) trends 
calculated over 1979-2005 for the past and 2006-2050 for 


the future (consistent with Son et al. [2010]) for the 
CCMVal-2 and CM1P5 simulations. The ozone trends 
amongst the different model groups (CMIP5 all, CHEM, 
NOCHEM, CCMVal) are comparable for the historical period. 
Error bars depict the ± 1 standard deviation across the model 
simulations and highlight the range of responses that 
result from interactively calculating ozone (CHEM and 
CCMVal-2) as opposed to prescribing it with a consistent 
database (NOCHEM). There is a clear relationship between 
the magnitude of the ozone depletion across the CHEM 
models and the near-global TLS trends (Figure 11a, 
represented here by the 100 hPa level; r= 0.60). However, this 
relationship breaks down in the future simulations where CO 2 
forcing becomes the dominant contributor to near-global 
lower stratospheric cooling, while the impact of ozone 
changes — the dominant contributor to past stratospheric 
cooling — reduces. 

[ 41 ] Antarctic polar stratospheric temperatures in spring 
(Figure 9b) are highly variable because they are affected 
not only by ozone changes but also by tropospheric wave 
forcing [Fogt et al., 2009]. RSS observations show a large 
decline of polar cap temperatures until 1999 followed by 
a 4 year period of higher temperatures in the smoothed 
time series. This rebound of temperature is likely episodic 
climate variability due to enhanced wave forcing from the 
troposphere [Langematz and Kunze, 2006], rather than a 
signal of ozone recovery, and is followed by considerably 
lower temperatures again after 2005 (not shown). The 
CMIP5 models simulate on average a slightly larger 
temperature decline (—1.29 ± 0.84 K/decade) than that 
observed (—0.79 ± 0.55 K/decade), but the observed trend 
is clearly within the spread of simulated trends. In the next 
section, we focus on the relationships between ozone, lower 
stratosphere temperatures, and the SH extratropical circula- 
tion in the model simulations. 

4.2. Ozone and the Southern Hemisphere 
Extratropical Circulation 

[ 42 ] Stratospheric ozone depletion in the Antarctic polar 
stratosphere also strongly impacts the SH extratropical circula- 
tion [Gillett and Thompson, 2003; Polvani et al., 2011b; 
Thompson and Solomon, 2002], Radiative cooling associated 
with ozone depletion leads to an increase in the meridional 
temperature gradient in the UTLS and a corresponding 
strengthening of the stratospheric polar vortex [ Thompson 
et al., 2011]. These stratospheric changes propagate 
down to the troposphere, reaching maximum amplitude in 
austral summer [Thompson and Solomon, 2002], While 
many previous modeling studies have found an association 
between ozone depletion and corresponding changes in 
the SH extratropical tropospheric circulation [Arblaster 
et al., 2011; Gillett and Thompson, 2003; McLandress 
et al., 2011; Polvani et al., 2011b], the mechanisms 
by which changes in the polar stratosphere impact the 
troposphere remain unclear. Various hypotheses have been 
suggested including changes to the low-level baroclinicity 
and planetary wave activity due to the strengthened polar 
vortex [Orretal., 2012; Simpson etal., 2012a] or via changes 
to synoptic wave activity which feeds back onto and results in 
meridional shifts in the mean flow [Chen and Held, 2007; 
Kushner and Polvani, 2004; Simpson et al., 2009]; see also 
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Figure 10. Trends in (a) annual-mean near-global (NG, — 82.5°S to 82.5°N) ozone at 50hPa, (b) annual 
mean near-global (82.5° S to 82.5°N) temperature at lOOhPa, (c) September-October-November-December 
(SOND) ozone at 50hPa over Antarctica (60°S-90°S), (d) October-November-December-January (ONDJ) 
temperature at lOOhPa over Antarctica (60°S-90°S), (e) DJF SH jet position at 850hPa, and (f) DIF upper 
tropospheric tropical (30°S-30°N) temperatures at 250 hPa.. Trends for the CM1P5 models are indicated 
by filled circles, CMIP5 NOCHEM models by triangles, CMIP5 CHEM models by squares, and CCMVal- 
2 by open circles, with different colors indicating the different experiments. The trends are calculated over 
1979-2005 for the past and over 2006-2050 for the future. 


Thompson et al. [201 1] and Gerber et al. [2012] for compre- 
hensive reviews and references therein. 

[ 43 ] During the late 20th century, ozone depletion acted in 
concert with increasing GHGs to impact the SH extratropical 
circulation and wider climate system [Son et al. 2010; 
Polvani et al., 2011a, 2011b; Kang et al. 2011]. However, 
ozone recovery will lead to increased temperatures in the 
Antarctic stratosphere in the future, thus having the opposite 
effect to ozone depletion. Previous studies have highlighted 
the importance of ozone recovery on either reversing or 
offsetting the impact of GHGs on the SH summertime 
atmospheric circulation projections [Arblaster et al., 2011; 
Perlwitz et al., 2008; Son et al., 2008; Son et al, 2010], with 
two recent studies [McLandress et al., 2012; Polvani et al., 
2011b] suggesting that the two forcings largely cancel out 
any impacts on the summertime atmospheric circulation over 


the coming half century under the SRES A1B GHG scenario. 
The exact balance between the impact of ozone recovery and 
increasing GHGs is both model and scenario dependent, with 
both the modeled rate of ozone recovery [Perlwitz, 2011] and 
their climate sensitivity [Arblaster et al., 2011] contributing to 
the magnitude of the atmospheric circulation response. 

[ 44 ] To provide some evaluation of the CM1P5 models, 
we first examine trends in December-January-February 
(DJF) zonal mean zonal winds over 1979-2005 for the 
CMIP5 multimodel mean compared to the ERA-Interim 
reanalysis data [Dee et al., 2011] in Figure 12. The CMIP5 
multimodel mean (Figure 12b) displays a robust poleward 
shift in the SH zonal mean zonal wind trends, consistent in 
magnitude with the reanalysis (Figure 12a). Figure 12 also 
depicts the climatological zonal mean zonal winds in thin 
black contours. An equatorward bias in the latitudinal 


5046 


EYRING ET AL.: OZONE AND ASSOCIATED CLIMATE IMPACTS 


(a) 



- 0.40 - 0.20 0.00 0.20 0.40 

1 00 hPa ANN NG T trend (K/dec) 



100 hPa ONDJ Temp (60-90S) trend 


(c) 



DJF jet position 850 hPa trend (deg/dec) 


(d) 


CMIP5 all =015, -014 
Chem CMIP5 = 0.20, -0.17 
noChemCMIP5 =0.60, -0.65* 


°? 

(h 


“ 0.2 - 


DJF jet position 850 hPa trend (deg/dec) 


A noohem CMIP5 □ chem CMIP5 O CCMVal-2 
H historical E3 REF-B2 1979-2005 rcp26 □ rcp45 B rcp60 ■ rcp85 D REF-B2 2006-2050 

Figure 11 . Trend relationship between (a) annual mean near-global mean ozone at 50 hPa and temperature 
at 100 hPa; (b) SOND ozone at 50hPa and ONDJ temperature at 100 hPa over Antarctica (60°S-90°S), 
(c) SOND ozone at 50hPa and DJF jet position at 850 hPa, and (d) DJF 250 hPa tropical (30°S-30°N) 
temperatures and DJF jet position at 850 hPa. The correlation between the two variables is calculated 
for all CMIP5 models and for the NOCFLEM and CHEM models for the past (first number) and all RCPs 
(second number). The correlation coefficients which are statistically significant at the 5% level 
(95% confidence) are indicated with an asterisk. The trends are calculated over 1979-2005 for the past and 
2006-2050 for the future. 


position of the SH extra tropical jet was found in the major- 
ity of CMIP3 models [ Kidston and Gerber, 2010], Swart 
and Fyfe [2012], Wilcox et al. [2012], and Ceppi et al. 
[2012] find a similar bias is evident in the CMIP5 models 
although somewhat reduced compared to CM1P3 models. 
This bias has implications for the response of the 
extratropical circulation to ozone and other externally 
forced changes such as GHGs [Bracegirdle et al., 2013; 
Kidston and Gerber, 2010; Wilcox et al., 2012], although 
the exact mechanism by which this impacts the circulation 
is a focus of current research [Simpson et al., 2012b]. 

[ 45 ] Future trends in zonal mean zonal winds calculated 
over the 2006-2050 period for the four RCP scenarios are 
shown in Figures 12c-12f for the CM1P5 multimodel 
mean. Although the radiative forcing differences over this 


near-term period are not large, it is clear that as the anthro- 
pogenic GHG forcing increases from RCP 2.6 to 8,5, so 
does the impact on the DJF SH extratropical jet, with 
robust poleward shifts under RCP 8.5. Under the high- 
mitigation scenario (RCP 2.6), which exhibits very little 
warming in the troposphere, ozone recovery appears to 
dominate the response with a weak equatorward shift of 
the SH jet, although the multimodel change is not signifi- 
cant. In RCP 4.5 and 6.0, the balance between increasing 
GHGs and ozone recovery largely offsets any changes to 
the SH circulation. Projected changes in the lower tropo- 
spheric winds (see Appendix Figure A3) are quite zonally 
symmetric in the SH extratropics and show weak equator- 
ward shifts in the SH midlatitude jet under RCP 2.6 and a 
robust poleward shift under RCP 8.5. 
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Figure 12. Long-term mean (thin black contour) and linear trend (color) of zonal mean DJF zonal winds 
for (a) ERA-Interim over 1979-2005; (b) CMIP5 (all models listed in Table 1 except those noted with b ); 
and (c, d, e, and f) same as Figure 12b but for the four RCPs over 2006-2050. Contour intervals of clima- 
tological wind are lOms” 1 starting from — 20ms _1 . Trends which are statistically significant at the 5% 
level (95% confidence) are stippled. 


[46] Figures 10 and 11 investigate the impact of 
stratospheric ozone on the SPl climate trends of the CFLEM, 
NOCHEM, CMlP5-all (both CHEM and NOCHEM), and 
CCMVal-2 models. CFIEM models that were excluded earlier 
due to ozone biases are also excluded here, and only 
NOCHEM models for which ozone forcing was explicitly 
available at the CMIP5 data archive are analyzed. Individual 
model trends are shown in Figure 1 1, with correlation coeffi- 
cients within each panel indicating the linearity of the relation- 
ships across the models for the historical and future periods. 

[ 47 ] The multimodel mean trends in springtime SH 
polar lower stratospheric ozone and temperatures 
(Figures 10c and lOd) indicate that, while the CMIP5 histori- 
cal and future changes are very similar in magnitude to the 
CCMVal-2 models in the multimodel mean, there is a large 


spread amongst the CHEM models, with somewhat weaker 
(and again, more consistent) trends for the NOCHEM models. 
Little variation in Antarctic ozone recovery or lower strato- 
spheric temperatures is found amongst the different RCPs 
(see section 3.2). Significant correlations are found between 
the polar stratospheric ozone and temperature changes when 
examining the individual model trends (Figure lib; r = 0.76 
and 0.69 across the CHEM-historical and CHEM-future, 
respectively). The CHEM models have both weaker and stron- 
ger temperature and ozone responses compared to the 
NOCHEM models in both the historical and future periods. 
Note that this small variation in ozone accounts for the weak 
NOCHEM intermodel correlations (r=— 0.21 and —0.26 for 
the past and future, respectively) here. The CCMVal-2 models 
are interspersed amongst the CM1P5 models. 
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[48] As noted above, these ozone-driven temperature 
changes lead to shifts in the position of the SH tropospheric 
jet (at 850 hPa), which are shown in Figure lOe for the 
multimodel mean and Figure lie for individual models. 
The position of the westerly jet is calculated as the latitude 
of maximum zonal mean zonal wind at 850 hPa as in Son 
et al. [2010]. The jet shifts poleward in most models for 
the historical period due to the combined effects of GHGs 
and ozone depletion. The correlation between the size of 
the ozone forcing and the size of the jet shift is significant 
at r = 0.80 for the CHEM-historical, indicating a strong 
sensitivity of the SH tropospheric circulation to ozone 
forcing. The NOCHEM-historical correlation is close to 
zero, highlighting the similarity of their ozone forcing. 

[ 49 ] The CMIP5 multimodel trends of the jet position 
under the four RCPs (Figure lOe) for the future projections 
until 2050 change from a weak equatorward shift for RCP 
2.6 to a weak poleward shift in RCP 8.5 (see also Figure 12) 
due to the increasing influence of the GHG forcing. The jet 
position trend is consistently more poleward in the NOCHEM 
than in the CHEM models, consistent with the weaker ozone 
changes in the NOCHEM models (Figure 10c). For the future 
projections, the picture is complicated by the opposing forces 
of ozone recovery and GHG warming on the jet position, with 
the relative impact of ozone recovery declining as the strength 
of the GHG forcing increases. Thus, the link between ozone 
forcing and jet position for future trends is reduced (compared 
to the past) with a correlation of r=0.35 for the CHEM 
models. The NOCHEM models again show little sensitivity 
to ozone (r= —0.14) in the future. 

[ 50 ] The impact of GHG warming on the jet position 
can be further identified by examining the variations of jet 
position with changes in tropical upper tropospheric temper- 
atures (Figures lOf and lid). Since ozone trends are very 
similar across the NOCHEM models, the jet position trend 
variations are primarily associated with variations in the 
response to GHGs, both in terms of the range of RCP forcing 
and the climate sensitivity of the models. The significant 
correlation across the NOCHEM models for the future 
period (r=— 0.68) indicates that the model response to 
GHGs also plays an important role in modifying future shifts 
in the SH tropospheric extratropical circulation, consistent 
with previous studies [Arblaster et al., 2011; Wilcox et al., 
2012]. Weaker correlations across the CHEM models 
(r=— 0.17 for the future) are likely due to the offsetting 
influence of ozone recovery, indicating that the rate of ozone 
recovery can be of similar importance to the response to 
GHGs, especially for the lower emission scenarios. Note that 
for the historical period the NOCHEM models show a 
positive, although insignificant, correlation. Other factors 
besides GHGs, such as decadal variability in ENSO, 
nonzonal impacts on the jet, or other forcing differences, 
could be driving the variations across the models for this 
period where GHG forcing is weaker than in the future. A 
weaker historical role for GHGs is also consistent with 
previous studies that argue for the dominance of ozone 
depletion on historical shifts in the jet [Polvani et al., 201 la] 

[ 51 ] Weak changes in the tropospheric jets under the RCP 
4.5 and 6.0 (Figure 12) are similar to that found for the 
CCMVal-2 simulations under the SRES A1B scenario [Son 
et al. [2010], Figure 4], The response in the stratosphere in 
the CMIP5 multimodel mean, however, is much weaker than 


in the CCMVal-2 mean. Although this is somewhat surpris- 
ing, given the two simulations have similar ozone forcing 
and the SRES A1B GHG forcing lies between RCP 
4.5 and 6.0, Figure 11 shows that the size and even sign of 
future trends in the position of the tropospheric jet to mid- 
century varies widely across models. 

[ 52 ] In summary, it is clear that ozone recovery plays an 
important role in the SH summertime extratropical circula- 
tion projections, with a reversal of the recent poleward shift 
in the jet possible under the high-mitigation scenario of RCP 
2.6. We find a strong sensitivity of the SH climate to ozone 
forcing, which is not captured by using models with a 
common prescribed ozone data set. This is an important 
complexity that needs to be considered in studies of SH 
climate change. We emphasize that the homogeneity of the 
ozone forcing data set across the CM1P5 model simulations 
assumed by Wilcox et al. [2012] for a limited number of 
models is not as clear-cut for the larger sample shown here. 

5. Summary 

[ 53 ] In this study, the representation of ozone and 
associated climate impacts in 46 model versions participating 
in the Coupled Model Intercomparison Project Phase 5 
(CM1P5) [ Taylor et al., 2012] project has been analyzed. We 
first evaluated the models' ability to reproduce ozone, zonal 
wind, and lower stratospheric temperature changes over the 
recent decades and then analyzed the future evolution of 
ozone from 2006 to 2100 under the four Representative 
Concentration Pathways (RCPs), again highlighting climate 
impacts associated with stratospheric ozone changes. 

[ 54 ] In CM1P3, around half of the models prescribed a 
constant ozone climatology, while the other half prescribed 
a time-varying ozone time series accounting for ozone 
depletion in the past and projected ozone recovery in the 
future [Son et al., 2008; Son et al., 2010], In contrast, 
all CMIP5 models considered time-varying ozone, either 
prescribed (28), interactively coupled (9), or calculated 
semi-offline (9), i.e., using a similar atmospheric model but 
prescribed SSTs and SICs and a detailed interactive chemis- 
try scheme. The main differences between the prescribed 
(NOCHEM) and interactive and semi-offline (CHEM) ozone 
is that the CHEM ozone projections are driven with the RCP 
GHG emission scenarios, whereas the ozone concentrations 
in the prescribed class are from a single GHG scenario 
different from the RCPs. For those using prescribed ozone, 
23 of the 28 applied the recommended IGAC/SPARC ozone 
database [Cionni et al., 2011], thus following the SRES A1B 
GHG scenario in the future which was used to drive the 
CCMVal-2 simulations. The remaining five used a different 
data set but still based on a single GHG scenario (see details 
in the Appendix). 

[ 55 ] Overall, ozone in the set of CMIP5 models analyzed 
here is reasonably well represented, but some large devia- 
tions exist for individual models, in particular, models with 
interactive chemistry. One goal of this paper was to docu- 
ment these differences for follow-up studies that require this 
information. For example, the recently published detection 
and attribution study by Santer et al. [2012] clearly illus- 
trates that it is very difficult to understand and interpret the 
causes of differences between modeled and observed atmo- 
spheric temperature trends without accurate information 
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regarding the forcings used in simulations of historical 
climate change. 

[56] Another goal was to test the sensitivity of ozone 
projections to GHG scenarios. The results found here 
confirm previous conclusions [Evring et al., 2010a], 
showing that future stratospheric ozone evolution is most 
sensitive to GHG scenarios in the Arctic (~40 DU difference 
among the four RCPs by 2100) and the northern midlati- 
tudes (-20DU difference among the four RCPs by 2100), 
with little sensitivity in the Antarctic ozone hole. For ozone 
in the troposphere, the projections broadly confirms previous 
studies in that, due to the projected reductions in NO x , VOC, 
and CO emissions, the future evolution of tropospheric 
column ozone is largely sensitive to the CH 4 level, together 
with an increased importance of the influx of ozone from the 
stratosphere [Kawase et al., 2011; Lamarque et al., 2011; 
Young et al., 2013b]. Other climate impacts on tropospheric 
ozone, such as increased loss through higher specific humid- 
ity, changes in lightning NO x emissions, and changes to 
reaction rates, could also be important but need to be isolated 
through sensitivity simulations and more detailed ozone 
budget diagnostic output. 

[ 57 ] Finally, we examined important climate impacts 
that are associated with stratospheric ozone changes. Our 
study confirms the findings of Santer et al. [2012] that the 
CMIP5 multimodel mean underestimates the annual near- 
global lower stratospheric temperature (TLS) trends with 
the NOCHEM models on average exhibiting a slightly larger 
trend than the CHEM models. Observed Antarctic polar 
temperature trends in spring are within the spread of CMIP5 
models. The analysis of temperature trends both globally 
and over the Antarctic illustrate the controlling factor of 
ozone changes on simulated past TLS temperature trends. 
In the NOCHEM models, in which ozone changes are 
prescribed with a very similar magnitude among the models, 
temperature trends show a smaller spread than CHEM 
model in which stratospheric ozone changes are simulated 
and affected by the interactions of stratospheric ozone 
chemistry and climate. 

[58] Significant impacts of stratospheric ozone depletion 
and recovery on the SH extratropical circulation in the 
CMIP5 models are also found, confirming previous studies 
with earlier generations of models [Son et al., 2008; Son 
et al., 2010], The influence of ozone is clearly seen in 
the CM1P5 multimodel mean, with a weak equatorward shift 
of the SH DJF midlatitude tropospheric jet under the high- 
mitigation scenario (RCP 2.6) when warming due to GHG 
is limited. However, a large range of circulation responses 
are found for the future, indicating a strong interplay between 
forcing from GHGs and ozone recovery. There is a signifi- 
cant relationship between the strength of future ozone and 
jet changes, highlighting the importance of future ozone 
changes on SH climate and the requirement to incorporate 
realistic ozone changes. 

[ 59 ] Overall, through the use of either a time-varying 
prescribed ozone data set or interactive chemistry, the 
CMIP5 models represent a notable improvement on many 
of their CMIP3 counterparts, especially as several studies 
have demonstrated the important climatic impacts of 
changes in the abundance of this species. However, with this 
added complexity likely comes increased model diversity 
and perhaps increased uncertainty in future projections 


[Stainforth et al., 2007], This underscores the need for 
well-thought out process-oriented diagnostics such as those 
developed by the Chemistry-Climate Model Validation 
(CCMVal) activity [Eyring et al., 2005; SPARC-CCMVal, 

2010] , as well as well-designed sensitivity experiments, both 
to characterize the model diversity and, together with obser- 
vations, identify key areas for improvement. It is therefore 
recommended to collect output for process-oriented ozone 
evaluation from climate models with interactive chemistry 
in future CMIP phases. 

Appendix A 

[ 60 ] In this Appendix, we provide descriptions of the 
chemistry schemes used in the CMIP5 models with interactive 
and semi-offline chemistry (section Al) and for the models 
that prescribed ozone (section A2). In addition, we show a 
figure similar to Figure 1 (total column ozone climatologies, 
1986 to 2005 mean) of the primary article for the individual 
CM1P5 models that are part of the CHEM group plus the 
models that use a prescribed ozone data set different from 
the IGAC/SPARC ozone database (Figure Al ). These include 
the following: 

[61 ] 1. The nine CMIP5 models with interactive chemistry: 
CESM 1 (FASTCHEM), CESMl(WACCM), CNRM-CM5, 
GFDL-CM3, GISS-E2-H-p2, GISS-E2-H-p3, GISS-E2-R-p2, 
GISS-E2-R-p3, and MIROC-ESM-CHEM (see description in 
section A 1.1) 

[ 62 ] 2. The two CM1P5 models with semi-offline chemistry: 
CCSM4 and IPSL-CM5A-LR (see description in section Al .2). 

[63] 3. The MIROC4h model which uses a prescribed 
ozone database different from the IGAC/SPARC ozone da- 
tabase (see description in section A2.1). The other model 
that used a different field was GISS-E2-H pi for which the 
files were not available on the Earth System Grid (ESG). 

[64] Figure A2 shows maps of tropospheric columns 
for selected decades in the past and future, and Figure A3 
shows trends in DJF zonal winds at 850 hPa calculated for 
2006-2050 for the four RCPs. 

Al. CM1P5 Models With Interactive 
or Semi-Offline Chemistry (CHEM) 

Al.l. Interactive Chemistry in the CM1P5 Simulation 

[65] In nine of the 46 model versions, ozone chemistry 
is calculated interactively in both the stratosphere and 
troposphere (CESM 1 (FASTCHEM), CESMl(WACCM), 
CNRM-CM5, GFDL-CM3, GISS-E2-H-p2, GISS-E2-H-p3, 
GISS-E2-R-p2, GISS-E2-R-p3, and MIROC-ESM-CHEM). 
The main characteristics of the chemistry schemes are briefly 
discussed for the different models below. 

Al.1.1. CESMl(FASTCHEM) 

[66] The Community Earth System Model (CESM) atmo- 
spheric chemistry-climate model is based on the Community 
Climate System Model version 4 (CCSM4) [Gent et al., 

2011] , configured with a fully interactive ocean, sea ice, 
land, and terrestrial carbon cycle. The CMIP5 simulations 
with interactive chemistry were performed at 1.875° x2.5° 
horizontal resolution, with 26 vertical layers extending to 
3.5 hPa. The implementation of the processes that involve 
interaction between chemical tracers and the physical atmo- 
sphere (such as wet deposition and convection) is described 
by Lamarque et al. [2012], The chemistry mechanism uses a 
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Figure Al. Same as Figure 1, but for individual CMIP5 models of the CHEM group plus MIROC4h 
which uses a prescribed ozone data set different from the IGAC/SPARC database. 


simplified version of the LLNL-IMPACT chemistry scheme 
\Rotman et al., 2004] to simulate the main free troposphere 
processes relevant to climate change [Cameron-Smith 
et al., 2006] but with addition of a single isoprene tracer 
and omitting the species that are only relevant to the strato- 
sphere. Stratospheric ozone is simulated instead using Linoz 


version 2 [Hsu and Prather, 2009], with an additional sim- 
ple loss to simulate polar ozone holes using the method of 
Cariolle et al. [1990], The whole chemical mechanism uses 
15 tracers. The 3-D methane distribution is specified using 
the output of simulations of an earlier model (CAM3.5) that 
used a more complete chemical mechanism but fixed SSTs/ 
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Figure A2. Same as Figure 4, but tropospheric ozone columns for the individual CMIP5 models com- 
pared to MLS/OMI. The numbers in the titles show the near-global (60°S-60°N) area-weighted mean. 


SICs [Lamarque et al., 2011]. The chemical emissions for 
the historical CM1P5 simulations followed Lamarque et al. 
[2010]. Due to a couple of mistakes, the CMIP5 simulations 
with this model were run with zero NO x from lightning and 
an incorrect rate constant for the CO + OF! reaction. The rate 
constant calculation inadvertently ignored the chemical 
activation channel due to a sign error in the original refer- 
ence [ Senosiain et al., 2003], which approximately halved 
the globally averaged rate constant. Fortunately, the model 
responded by producing excessive CO concentrations that 
almost exactly compensated for the incorrect rate constant, 
such that the reaction rate (given by k[CO][OFl]) was 
relatively unchanged, and that error had little effect on the 
rest of the simulation. Specifically, the rate constant error 
only reduced globally averaged ozone by around 1 pbb and 
reduced ozone in highly polluted regions by 2-4 ppb. 
Al.1.2. CESMl(WACCM) 

[67] CESM 1 (WACCM) is a configuration of the NCAR 
Community Earth System Model (CESM1) with the 
Whole Atmosphere Community Climate Model (WACCM) 
version 4 as its atmospheric component \Calvo et al., 2012], 


WACCM has 66 vertical levels from the surface to 

4.5 x 10 6 hPa and runs at horizontal resolution of 1.9° 

2.5 coupled to the interactive ocean and sea-ice model 
components [Danabasoglu et al., 2012; Holland et al., 
2012], The capability to run WACCM as the atmospheric 
component in coupled simulations is one of the significant 
improvements over the immediate predecessor to CESM, 
the CCSM4 [Gent et al., 2011]. The CMIP5 simulations 
used WACCM version 4, which is a configuration of the 
Community Atmospheric Model, version 4 (CAM4), and 
includes all of the physical parameterizations of that model 
[Neale et al., 2013]. WACCM4 is a fully interactive model, 
where the radiatively active gases (C0 2 , ECO, N 2 0, CFFt, 
CFCs, NO, O 3 ) affect heating and cooling rates and there- 
fore dynamics. The model uses a finite volume dynamical 
core [Lin et al., 2004], which explicitly calculates the mass 
fluxes in and out of a given model grid cell, thus ensuring 
mass conservation. In terms of chemistry, WACCM4 
incorporates a detailed neutral mechanism for the middle 
atmosphere [Garcia et al., 2007]. The chemistry module is 
based on Model for Ozone And Related chemical Tracers 
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Figure A3. Zonal mean zonal wind DJF trends at 850 hPa for the CM1P5 multimodel mean over 2006-2050 
for (a) RCP 2.6, (b) RCP 4.5, (c) RCP 6.0, and (d) RCP 8.5. Model agreement is assessed through testing for 
significance in each model: stippling indicating most models show significant change and agree on the sign. 


(MOZART) version 3 [Kinnison et al., 2007] and includes 
O x , NO x , HO x , C10 x , and BrO x chemical families, along 
with CH 4 and its degradation products (a total of 59 species 
and 217 gas-phase chemical reactions). This version of 
WACOM does not include a detailed representation of 
tropospheric chemistry. Rate constants are based on Sander 
etal. [2006]. In addition, there are 1 7 heterogeneous reactions 
on three aerosol types: nitric acid trihydrate, supercooled 
ternary solution, and water ice. The reaction scheme, rates, 
and method of solution are identical to those evaluated for 
the SPARC-CCMVal [2010] report. 

Al.1.3. CNRM-CM5 

[68] The atmospheric component of the CNRM-CM5.1 
model [Voldoire et al., 2012] is the global spectral 
ARPEGE-Climat atmospheric model. All CMIP5 simula- 
tions were performed at 1.4° x 1.4° horizontal resolution, 
with 31 vertical layers extending to lOhPa, including an 
interactive linear ozone chemistry. Within ARPEGE-Climat, 
the ozone mixing ratio is treated as a prognostic variable 
with photochemical production and loss rates computed by a 
two-dimensional zonal chemistry model (MOBIDIC, Cariolle 
and Teyssedre [2007]) and is used for the radiative transfer 
calculation. In this scheme, the net photochemical production 
in the ozone continuity equation is expanded into a first- 
order Taylor series around the local value of the ozone mixing 
ratio, the temperature, and the overhead ozone column. An 
additional ozone destruction term is used to parameterize the 
heterogeneous chemistry as a function of the equivalent 
chlorine content prescribed for the actual year. The linearized 
scheme is based on a full chemistry scheme with 56 constitu- 
ents, 175 chemical reactions, and 51 photo-reactions, 
implemented in a two-dimensional stratospheric model. 
The results of a century integration with this two-dimensional 


model have been used for the linearization of the photochemis- 
try near the equilibrium. The error between the full and the lin- 
earized scheme is smaller than 10% in the two-dimensional 
model. The ozone mixing ratio remains a prognostic variable 
in the troposphere, but the parameterization is mainly devel- 
oped for use in the stratosphere (i.e., some tropospheric chem- 
ical reactions are not taken into account). 

Al.1.4. GFDL-CM3 

[69] In the atmospheric component of the GFDL-CM3 
model (called AM3) [Donner et al., 2011], the chemistry 
models of Horowitz et al. [2003] for the troposphere and 
Austin and Wilson [2006] for the stratosphere are merged. 
The CM1P5 simulations with interactive chemistry were 
performed on a cubed-sphere grid [Putman and Lin, 2007] 
at C48 horizontal resolution (with grid cell sizes varying 
from 163 km to 231km), with 48 vertical layers extending 
to 0.01 hPa. Vertical grid spacing is approximately 1.5 km 
near the tropopause increasing to about 3 km near the 
stratopause. AM3 simulates the atmospheric concentrations 
of 97 chemical species (81 gases and 16 aerosols), which 
undergo 172 gas-phase reactions, 41 photolytic reactions, 
and 11 heterogeneous reactions. Stratospheric chemistry 
includes O x , NO x , HO x , C10 x , and BrO x chemical species 
[Austin and Wilson, 2006], The source of reactive chlorine 
and bromine are parameterized as a function of tropospheric 
concentrations of source gases (lagged by the stratospheric 
“age of air”) [Austin and Wilson, 2006], The parameteriza- 
tion uses observed source gas distributions to estimate 
the fractional rate at which source gases entering the strato- 
sphere are photolyzed and converted into reactive halogen 
species along their transport path through the stratosphere. 
Heterogeneous reactions are included on ice and nitric acid 
trihydrate polar stratospheric clouds (PSCs) and in liquid 
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ternary solution aerosols [Carslaw et al., 1995; Hanson 
and Mauersberger , 1988]. Tropospheric chemistry includes 
O x -NO x -HO x -CO-CH 4 -NMHC hydrocarbon chemistry 
(including alkanes up to a lumped 4+ carbon species, 
alkenes up to propene, and isoprene and terpenes). Clear-sky 
photolysis frequencies are calculated using a multivariate 
interpolation table derived from the tropospheric ultraviolet- 
visible radiation model [Madronich et al., 1998] suitable for 
calculating tropospheric and stratospheric photolysis rates, 
with an adjustment applied for the effects of large-scale clouds 
[. Brasseur et al., 1998]. Kinetic reaction rates and mass 
accommodation coefficients and reaction probabilities are 
taken from Sander et al. [2006]. 

Al.1.5. GISS-E2 

[70] The new version of the GISS climate model used 
for the CMIP5 simulations is called ModelE2. It is similar 
to ModelE [ Schmidt et al., 2006] but has numerous improve- 
ments to the physics, including the ability to represent 
multiple downdrafts and updrafts in convective systems 
( Shindell et al., in prep., 2012). The CMIP5 simulations with 
interactive chemistry were performed at 2° x 2.5° horizontal 
resolution, with 40 vertical layers extending to 0. 1 hPa. 
The atmosphere was coupled to a full dynamic ocean, either 
the “Russell” ocean [ Russell et al., 1995] (GISS-E2-R) 
or the Hybrid Coordinate Ocean Model (HYCOM) ocean 
[Sun and Bleck, 2006] (GISS-E2-H). The indices p in Table 1 
refer to non-interactive aerosols/chemistry (p=l), interac- 
tive aerosols/chemistry (p = 2), and same asp = 2 but with 
an a priori calculation of the aerosol indirect effect instead 
of a parameterized scheme (p = 3). The GISS chemical 
mechanism is fully embedded within ModelE2, so that 
chemical constituents are treated consistently with the 
physics of other parts of the model, such as surface fluxes 
of fundamental physical quantities (e.g., heat) and with 
transport of momentum and other constituents such as water 
vapor. To represent stratospheric chemistry, the model 
includes chlorine- and bromine-containing compounds 
and CFC and N 2 0 source gases. PSC formation in the 
stratosphere is dependent upon water vapor, temperature, 
and HNO3 [Hanson and Mauersberger, 1988], and the 
heterogeneous hydrolysis of N2O5 on sulfate follows Kane 
et al. [2001] and Hallquist et al. [2003], Chemistry and 
aerosols are fully integrated, so that these components 
interact with each other and with the physics of the climate 
model. The model contains 51 species for gas-phase 
chemistry interacting via 156 reactions. Tropospheric 
chemistry includes basic NO x -HO x -O x -CO-CH 4 chemistry 
as well as Peroxyacetyl nitrate (PAN) and the hydrocarbons 
isoprene, alkyl nitrates, aldehydes, alkenes, and paraffins. 
Acetone chemistry [Houweling et al., 1998] and a newly 
identified reaction channel of OH + N0 2 to fonn HNO3 
[Butkovskaya et al., 2007], the latter with a reaction rate 
coefficient that depends on pressure and temperature, has 
been added. 

Al.1.6. MIROC-ESM-CHEM 

[71] MIROC-ESM-CHEM is a version of the Earth 
system model MIROC-ESM [Watanabe et al., 2011] 
coupled with atmospheric chemistry. The atmospheric 
chemistry component is based on that used in the chemistry 
model CHASER [Sudo and Akimoto, 2007; Sudo et al., 
2002], which considers the detailed photochemistry in the 
troposphere, wet and dry deposition of chemical species, 


and emissions of primary chemical species from a variety 
of sources. For Model for Interdisciplinary Research on 
Climate (MIROC)-ESM-CHEM, the chemistry of CHASER 
has been extended to include the stratosphere by incorporat- 
ing halogen chemistry and related processes based on the 
schemes used in the Climate System Research/National 
Institute for Environment Studies (CCSR/NIES) stratospheric 
chemistry model [Akiyoshi et al., 2004]. In the configuration 
used for the CMIP5 simulations, the model calculates the 
concentrations of 92 chemical species with 262 chemical 
reactions at approximately 2.8° by 2.8° (T42) horizontal 
resolution with 80 vertical layers up to 0.003 hPa. The 
model considers the O x -NO x -HO x -CH 4 -CO chemistry 
together with the oxidation of several non-methane VOCs 
(ethane, propane, propene, butane, acetone, methanol, 
isoprene, and terpenes) to represent the ozone chemistry in 
the troposphere. The model has detailed stratospheric 
chemistry, including reactions of chlorine- and bromine- 
containing compounds, CFCs, HFCs, OCS, and N 2 0, as well 
as the formation of PSCs and associated heterogeneous reac- 
tions on their surfaces. Photolysis rates (J- values) are calcu- 
lated online using temperature and actinic fluxes computed 
in the radiation component [Sekiguchi and Nakajima, 2008] 
considering absorption and scattering by gases, aerosols, 
and clouds as well as the effect of surface albedo. The influ- 
ences of short-wave radiative forcing (RF) associated with 
the solar cycle, volcanic eruptions, and subsequent changes 
in stratospheric ozone are also taken into account for the 
calculation of the photolysis rate. 

A1.2. Semi-Offline Chemistry in the CMIP5 Simulation 

[72] In nine of the 46 model versions, ozone chemistry 
was calculated semi-offline in both the stratosphere and 
troposphere (BNU-ESM, CCSM4, CESMl(BGC), CESM1 
(CAM5), IPSL-CM5A-LR, IPSL-CM5A-MR, IPSL-CM5B- 
LR, NorESMl-M, NorESMl-ME). These semi-offline ozone 
data sets are calculated in two different ways, which are 
described in the two subsections below. In the ozone figures 
that show individual models, only one representative of each 
of the two data sets is included (CCSM4, IPSL-CM5A-LR). 
Similarly, only CCSM4 is included in the CHEM ozone 
multimodel mean together with the models with interactive 
chemistry (section A 1.2), while 1PSL-CM5A is excluded 
because it uses a simplified scaling to estimate GHG depen- 
dence of stratospheric ozone. 

Al.2.1. Community Atmosphere Model 
With Interactive Chemistry 

[73] In a subset of the model simulations (BNU-ESM, 
CCSM4, CESMl(BGC), CESM1(CAM5), NorESMl-M, 
and NorESMl-ME, entry P 2 in Table 1), ozone is calculated 
semi-offline using the global 3-D Community Atmosphere 
Model 3.5 (CAM3.5) version [Gent et al., 2011] with 
interactive chemistry to calculate distributions of gases and 
aerosols in the troposphere and the lower to midstratosphere 
[Lamarque et al., 2011; Lamarque et al., 2012; Lamarque 
et al., 2010]. CAM3.5 also participated in CCMVal-2 and 
was evaluated as part of the SPARC-CCMVal [2010] report. 
In order to limit computational cost, this model only solves 
for the atmospheric and land portions of the climate system, 
using precomputed SSTs and SICs as boundary conditions. 
The model configuration consists of a horizontal resolution 
of 1.9° x 2.5° and 26 hybrid levels, from the surface to 
3.5 hPa with a 30 min timestep. The chemical mechanism 
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used is formulated to provide an accurate representation of 
both tropospheric and stratospheric chemistry. Specifically, 
a representation of stratospheric chemistry (including 
polar ozone loss associated with stratospheric clouds) from 
MOZART version 3 [Kinnison et ah, 2007] is included. 
The tropospheric chemistry mechanism has a simplified 
representation of non-methane hydrocarbon chemistry in 
addition to standard methane chemistry, extended from 
Houweling et al. [1998] with the inclusion of isoprene and 
terpene oxidation. A monthly ozone climatology is created 
by averaging monthly mean ozone fields from the transient 
simulations for calendar months over one decade. These 
climatological fields are read as 3-D fields used in the radia- 
tive transfer calculation in the Community Climate System 
Model version 4 (CCSM4; Meehl et al. [2012]) and in the 
models listed above. The ozone time series was provided 
as decadal-average monthly mean concentrations, linearly 
interpolated to yearly resolution within the coupled models. 
These model simulations are in the class of interactive 
and semi-offline chemistry (CHEM) in this paper since 
stratospheric ozone follows the RCPs in the future unlike 
the IGAC/SPARC database that derives stratospheric ozone 
concentrations based on the SRES A1B GHG scenario. 
Al.2.2. IPSL 

[ 74 ] In the three IPSL coupled ocean-atmosphere climate 
model versions (IPSL-CM5A-LR, IPSL-CM5A-MR, and 
1PSL-CM5B-LR, [Dufresne et al., 2013]), ozone fields are 
calculated semi-offline by a tropospheric and a stratospheric 
version of the atmospheric general circulation model LMDz 
with interactive chemistry, following the four RCPs in the 
future. Tropospheric ozone is calculated using the tropo- 
spheric chemistry module INCA (Interaction with Chemistry 
and Aerosol) within a version of LMDz that has 19 vertical 
levels, with most of them in the troposphere [ Hauglustaine 
et al., 2004]. Stratospheric ozone is calculated using the 
stratospheric chemistry module REPROBUS (Reactive 
Processes Ruling the Ozone Budget in the Stratosphere), 
using a stratospheric version of LMDz that has 50 vertical 
levels, with a resolution varying in the stratosphere from 
about 1 km around the tropopause region to about 3 km at 
50 km \Jourdain et al., 2008], 

[ 75 ] In the case of LMDz-REPROBUS, only two simu- 
lations for the future were performed: a CCMVal REF- 
B2 simulation that follows the SRES A1B GHG scenario 
and a CCMVal SCN-B2c simulation, where GHGs were 
fixed at 1960 levels throughout the simulation [Eyring 
et al., 2010b]. For this reason, a rather crude approxima- 
tion was used based on previous calculations carried out 
with CAM3.5, the only CCM that ran all the RCP calcu- 
lations in Eyring et al. [2010a], For CAM3.5, correlation 
plots for different vertical and latitudinal regions of 
the stratosphere indicate a somewhat linear dependency 
of stratospheric ozone change on CO 2 level change 
over the range of RCP scenarios (see Figure 3 of Eyring 
et al. [2010a]). Assuming a similar linear dependency in 
the case of LMDz-REPROBUS, stratospheric RCP 2.6, 
RCP 4.5 and RCP 8.5 ozone fields were generated by 
interpolating or extrapolating linearly from the REF-B2 and 
SCN-B2c stratospheric ozone time series. Ozone fields 
are reconstructed using locally derived time-varying 
weighting coefficients that are assumed to be proportional to 
the CO 2 mixing ratio. 


[76] Since atmospheric measurements and model 
simulations constrained by emission inventories indicate that 
the increase in Ozone Depleting Substances (ODSs; mostly 
CFCs) concentrations was very marginal before 1960 
[WMO, 2011], trends in stratospheric ozone are neglected 
before 1960. Raw monthly 3-D fields of both INCA and 
REPROBUS ozone time series are then zonally averaged 
and temporally smoothed using an 1 1 year running mean. 
Finally, the tropospheric and stratospheric zonal mean ozone 
time series extracted from the LMDz-INCA and LMDz- 
REPROBUS simulations are merged at the chemical 
tropopause (defined as the 150 nmol/mol ozone contour). 
INCA ozone fields are slightly stretched or compressed in 
order to have the INCA tropopause altitude matching the 
REPROBUS tropopause altitude. The thickness of the 
transition region is taken to be about 3 km, extending from 
2 km above the tropopause to 1 km below the tropopause. 
The resulting merged zonal mean ozone fields contain 
INCA fields below the tropopause region and REPROBUS 
fields above, together with a relatively smooth transition 
from one field to another within the tropopause region. 
The entire merging procedure is described in detail by Szopa 
et al. [2012], 

A2. CMIP5 Models With Prescribed 
Chemistry (NOCHEM) 

A2.1. IGAC/SPARC Ozone Database 

[ 77 ] In 23 model versions, the original or a modified 
version of the IGAC/SPARC ozone database (see section 1) 
was used. To represent this class of models in the figures, only 
the IGAC/SPARC ozone database is shown in the figures. 
The original IGAC/SPARC ozone database was recommended 
on the CMIP5 website as forcings for models without interac- 
tive chemistry. It was used in 13 of the 28 model versions with 
prescribed ozone (ACCESS1-0, ACCESS1-3, BCC-CSM1.1, 
BCC-CSM1.1-M, CSIRO-Mk3-6-0, EC-EARTH, FGOALS- 
g2, FGOALS-s2, FIO-ESM, GFDL-ESM2G, GFDL- 
ESM2M, INM-CM4, and MRI-CGCM3). Slightly modified 
versions were used in 10 of the 28 model versions with pre- 
scribed ozone (CanAM4, CanCM4, CanESM2, CMCC-CM, 
HadCM3, HadGEM2-A, HadGEM2-CC, HadGEM2-ES, 
MPI-ESM-LR, and MPI-ESM-P). 

[78] Since the solar cycle was not included in the future 
part of the IGAC/SPARC ozone data set, a modified version 
that included the effects of the solar cycle on future ozone 
concentrations was created by the Met Office [Jones et al., 
2011] and used by 7 of the 31 model versions with prescribed 
chemistry (CMCC-CM, HadCM3, HadGEM2-A, HadGEM2- 
CC, HadGEM2-ES, MPI-ESM-LR, and MPI-ESM-P). The 
solar regression index was used to build a future time series 
consistent with a repeating solar irradiance modeled as a 
sinusoid with a period of 1 1 years, with mean and maximum- 
minimum values corresponding to solar cycle 23 normalized 
against the 180.5 run time series used in the historical ozone 
data set. Future ozone was represented as: 03(t) = a*SOL + 
b*EESC + seasonal cycle + residuals. For consistency, the 
indices of the 1 1 year solar cycle (SOL) and total EESC were 
identical to those used to prepare the original data set. In the 
HadGEM2-ES model [Collins et al., 2011] tropospheric chem- 
istry was fully interactive, taking account of surface emissions 
of CH,, CO, NO x , HCHO, C 2 H 6 , C 3 H 8 , Me 2 CO, and MeCHO. 
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[ 79 ] CanAM4, CanCM4, and CanESM2 again used a 
slightly modified version of the IGAC/SPARC data set. Over 
the historical period (1850-2009), it is identical to the original 
database, except that it is zonally averaged. From 2010 
onwards, ozone in the troposphere is also identical to the 
original database for all RCPs, although zonally averaged, 
while ozone in the stratosphere has been modified. During this 
future period, stratospheric ozone was derived by combining 
two separate terms derived from multiple linear regression 
(MLR). The first term represents EESC from the historical 
period, extended into the future using predictions of EESC 
from the suite of the CCMVal-2 models. The CCMVal-2 time 
series of future ozone was independently fitted to an EESC and 
a linear trend using MLR. The linear trend term from the 
models was the second tenn used to construct the future ozone, 
combined with the EESC term described above. Construction 
of the ozone time series in this manner helped ensure the 
response of ozone to EESC in the historical period was equal 
to the response of ozone to EESC in the future, as halogen 
concentrations decrease. 

A2.2. Other Ozone Data Sets 
[so] In the remaining 5 of the 28 model versions with 
prescribed ozone (GISS-E2-H pi, GISS-E2-R pi, M1ROC- 
ESM, MIROC4h and M1ROC5), data sets different from 
the Cionni el al. [2011] database were used. In the ozone 
figures that show individual models, only one representative 
of the two different data sets are shown with dotted lines 
(GISS-E2-H pi, and MIROC4h). 

A2.2.1. MIROC-ESM, MIROC4h and MIROC5 
[si] MIROC-ESM and MIROC5 prescribed past and future 
ozone concentrations using the simulations of Kawase et al. 
[2011], whereas MIROC4h used the Sudo et al. [2003] ozone 
data set for past simulations and Kawase et al. [2011] for 
future simulations. Kawase et al. [2011] simulated the long- 
term evolution of tropospheric ozone using the chemistry- 
climate model CHASER [Sudo et al., 2002], driven by 
historical and future RCP scenarios. Stratospheric ozone was 
based on the late 20th century climatology of the SUNYA data 
set [Wang et al., 1995], scaled by the observed stratospheric 
ozone trend from Randel and Wu [2007] for the past and by 
the projected time evolution of EESC in the future simulation 
(following the halogen scenario Al of WMO [1999]). Strato- 
spheric ozone in this data set does not reflect future changes 
in ozone transport that could be induced by increases in 
GHGs, as is projected in the IGAC/SPARC ozone database. 
As a consequence, this data set has a straightforward recovery 
of stratospheric ozone in the future toward chemically 
unperturbed past state (e.g., ozone is similar between the late 
1 9th and late 2 1 st century, which are two periods with negligi- 
ble effects of ODSs). There is neither super-recovery nor 
ozone reduction in the tropical lower stratosphere in the future, 
which is seen in the IGAC/SPARC database. Effects of the 
1 1 year solar cycle and volcanic aerosols are not considered 
in this ozone data set. 

A2.2.2. GISS-E2-H pi, GISS-E2-R pi 
[ 82 ] In the GISS-E2 non-interactive chemistry simulations 
{p = I ), offline ozone fields were used from the prior genera- 
tion (CMIP3) GISS simulations [Hansen et al., 2007]. 
Tropospheric ozone changes were from an older chemistry- 
climate model [Shindell et al., 2003] driven by prescribed 
changes in precursor emissions, ocean boundary conditions 
(SSTs and SICs), and GHG abundances. Stratospheric ozone 


changes from 1979 to 1997 were from the observational 
analyses of Randel and Wu [1999]. Tropospheric ozone 
was kept constant from 1990 on, while stratospheric ozone 
was constant from 1997 to 2003 and then linearly returned 
to its 1979 value between 2004 and 2050. 
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